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Non-intersecting squared Bessel paths and multiple orthogonal 
polynomials for modified Bessel weights 

O . 

O . A.B.J. Kuijlaars, A. Martmez-Finkelshtein, and F. Wielonsky 

O 

Q 

Q>^ ■ We study a model of n non-intersecting squared Bessel processes in the confluent case: 

all paths start at time i = at the same positive value x = a, remain positive, and are 
conditioned to end at time t = T at x ^ 0. In the limit n — > oo, after appropriate rescaling, 
the paths fill out a region in the te-plane that we describe explicitly. In particular, the 
paths initially stay away from the hard edge at a; = 0, but at a certain critical time t* the 
(— I , smallest paths hit the hard edge and from then on are stuck to it. For t ^ t* we obtain the 

usual scaling limits from random matrix theory, namely the sine, Airy, and Bessel kernels. 
A key fact is that the positions of the paths at any time t constitute a multiple orthogonal 
polynomial ensemble, corresponding to a system of two modified Bessel-type weights. As 
a consequence, there is a 3 x 3 matrix valued Riemann-Hilbcrt problem characterizing this 
model, that we analyze in the large n limit using the Deift-Zhou steepest descent method. 
^ I There are some novel ingredients in the Riemann-Hilbert analysis that are of independent 

■ interest. 

m 
(N 



1 Introduction 



I Determinantal point processes are of considerable current interest in probability theory and 

■ mathematical physics, since they arise naturally in random matrix theory, non-intersecting 

paths, certain combinatorial and stochastic growth models and representation theory of large 
^ I groups, see e.g. Deift [22], Johansson [31], Katori and Tanemura [37], Borodin and Olshanski 

^ ■ [11], and many other papers cited therein. See also the surveys of Soshnikov [i^, Konig [38] . 

Hough et al. [30], and Johansson [32]. 

A determinantal point process is characterized by a correlation kernel K such that for 
every m the m-point correlation function (or joint intensities) takes the determinantal form 

det [K{xj,Xk)]j^k^^^^^^^^ 

We will only consider determinantal point processes on M. 

As pointed out by Borodin [9] certain determinantal point processes arise as biorthogonal 
ensembles, i.e., joint probability density functions on of the form 

V{xi,...,Xn) = -^det[fj{xk)]j^k=i,...,n(iet[gj{xk)]j,k=i,...,n (1-1) 
for certain given functions /i ,...,/„ , and gi, . . . ,gn- The correlation kernel is then given by 

n 

K{x,y) = Y,<Pjix)i^jiy) (1-2) 



1 



where (pj, tpj, j = 1, . . . ,n are such that 

span{</)i, ...,(/)„} = span{/i, . . . , spanlV'i, • • • , ^An} = spanl^i, . . . ,5„} 

and they have the biorthogonahty property 

/ (f>j{x)'ilJk{x)dx = 6j^k- 

JR 

The joint probabihty distribution function for the eigenvalues of unitary invariant ensem- 
bles of random Hermitian matrices {l/Zn)e~'^^^^^dM has the form where 

fj{x) = g,{x) = x^-ie-^^(-), j = 1, 2, . . . , n. (1.3) 

Orthogonalizing the functions (|1.3p leads to 

cl)j{x) = il)j{x) =pj_i(x)e"5^(^), j = 1,2,... ,n, 

where Pj-i is the orthonormal polynomial of degree j — 1 with respect to the weight e"^*-^) 
on M. The kernel ()1.2p is then the orthogonal polynomial kernel, also called the Christoffel- 
Darboux kernel because of the Christoffel-Darboux formula for orthogonal polynomials, and 
the ensemble is called an orthogonal polynomial ensemble |38j . 

Other examples for biorthogonal ensembles arise in the context of non-intersecting paths 
as follows. Consider a one-dimensional diffusion process X[t) (i.e., a strong Markov process 
on R with continuous sample paths) with transition probability functions pt{x,y), t > 0, 
X, y S M. Take n independent copies Xj{t), j = 1, . . . ,n, conditioned so that 

• -'^^^■(O) = Oj, Xj(T) = bj, where T > 0, and ai < 02 < • • • < a„, bi < b2 < ■ ■ ■ < bn are 
given values, 

• the paths do not intersect for < t < T. 

It then follows from a remarkable theorem of Karlin and McGregor [33] that the positions of 
the paths at any given time t G (0, T) have the joint probability density (II. ip with functions 

fj{x) =pt{aj,x), gj{x) = pT~t{x,bj), j = l,...,n. 

[Properly speaking the joint probability density function is first defined for ordered n-tuples 
xi < X2 < ■ ■ ■ < Xn only. It is extended in a symmetric way to all of M".] 

An important feature of determinantal point processes is that they seem to have universal 
limits. By now, this is well-established for the eigenvalue distributions of unitary random 
matrix ensembles. Indeed if is the eigenvalue correlation kernel for the random matrix 
ensemble (note the n-dependence of the potential) 

then we have under mild assumptions on V that 

lim —Kn.{x,x) =: p{x) 
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exists. In addition if V is real analytic, and if x* is in the bulk of the spectrum (i.e., p{x*) > 0), 
then (see [21]) 

^ f * ^ * y \ smTr(x — y) . ^. 

n^oo np{x*) "\ np{x*)^ np{x*) J Tr{x — y) 

Universality of local eigenvalue statistics is expressed by (jl.4p in the sense that the sine kernel 
arises as the limit regardless of V and x* . The universality (|1.4p is extended in many ways 
and (as its name suggests) under very mild assumptions (see the recent works |43[ ) . 

The limit (jl.4p does not hold at special points x* of the spectrum where p{x*) = 0. 
However it turns out that Kn has scaling limits at such special points that are determined by 
the macroscopic nature of x* , and in that sense they are again universal (see e.g. [14^ 115^ [T6\ 

mm)- 

It is reasonable to expect that such universal limit results hold generically for non- 
intersecting paths as well, although results are more sparse. For recent progress related 
to discrete random walks, random tilings and random matrices with external source see 

It is the aim of this paper to study a model of n non-intersecting squared Bessel processes 
in the limit n oo. Recall that if {X(t) : t > 0} is a d-dimensional Brownian motion, then 
the diffusion process 



R{t) = ||x(t)||2 = VMt? + ■■■ + Xd{t)^ t>o, 

is the Bessel process with parameter a = f — 1, while i?^(t) is the squared Bessel process 
usually denoted by BESQ'^ (see e.g. [Ml Ch. 7], [39] )• These are an important family of 
diffusion processes which have applications in finance and other areas. The well known Cox- 
Ingersoll-Ross (CIR) model in finance describing the short term evolution of interest rates 
or different models of the growth optimal portfolio (GOP) represent important examples of 
squared Bessel processes [29lll8]. The Bessel process R{t) for d = 1 reduces to the Brownian 
motion reflected at the origin, while for d = 3 it is connected with the Brownian motion 
absorbed at the origin [36l [37] . 

A system of n particles performing BESQ*^ conditioned never to collide with each other 
and conditioned to start and end at the origin, can be realized as a process of eigenvalues of 
a hermitian matrix-valued diffusion process, known as the chiral or Laguerre ensemble, see 
e.g. [271 ESI (Ml [M] and below. In this paper we consider the case where all particles start at 
the same positive value a > and end at 0. Of particular interest here is the interaction of 
the non- intersecting paths with the hard edge at 0. Due to the nature of the squared Bessel 
process, the paths starting at a positive value remain positive, but they are conditioned to 
end at time T at 0. After appropriate rescaling we will see that in the limit n — > oo the 
paths fill out a region in the tx-plane. The paths start at t = and initially stay away from 
the hard edge at x = 0. At a certain critical time the smallest paths hit the hard edge and 
from then on are stuck to it. The phase transition at the critical time is a new feature of the 
present model. It is a new soft-to-hard edge transition. 

We are able to analyze the model in great detail since in the confluent case the biorthogonal 
ensemble reduces to a multiple orthogonal polynomial ensemble, as we will show in Subsection 
[2] below. The correlation kernel for the multiple orthogonal polynomial ensemble is expressed 
via a 3 X 3 matrix- valued Riemann-Hilbert (RH) problem [6l [20] . 

We analyze the RH problem in the large n limit using the Deift-Zhou steepest descent 
method for RH problems [26]. There are some novel ingredients in our analysis which we 
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feel are of independent interest. First of all, there is a first preliminary transformation which 
makes use of the explicit structure of the RH jump matrix. It contains the modified Bessel 
functions Iq, and la+i and we use the explicit properties of Bessel functions. A result of the 
first transformation is that a jump is created on the negative real axis, see Section [3l 

The multiple orthogonal polynomials for modified Bessel functions were studied before 
by Coussement and Van Assche \18\ [T9] . We use their results to make an ansatz about an 
underlying Riemann surface that allows us to define the second transformation in the steepest 
descent analysis in Section [H The use of the Riemann surface is similar to what is done in 
[71 142] . In the appendix we mention an alternative approach via equilibrium measures and 
associated functions. The further steps in the RH analysis follow the general scheme laid 
out by Deift et al. [MIES] in the context of orthogonal polynomials. An important feature of 
the present situation is that there is an unbounded cut along the negative real axis and we 
have to deal with this technical issue in the construction of the global parametrix in Section 
m The construction of the local parametrices at the hard edge also presents a new technical 
issue, see Section [H 

The main results of the paper are stated in the next section. 



2 Statement of results 

2.1 Squared Bessel processes 

The transition probability density of a squared Bessel process with parameter a > — 1 is given 
by (see [HI [39]) 

Pf(x,y) = ^ (4^) , x,y> 0, (2.1) 

where la denotes the modified Bessel function of the first kind of order a, 

see [H Section 9.6] for the main properties of the modified Bessel functions. If d = 2(a + 1) is 
an integer, then the squared Bessel process can be seen as the square of the distance to the 
origin of a d-dimensional standard Brownian motion. 

If the starting points aj and the endpoints hj are all different, then (as explained in the 
introduction) the positions of the paths at a fixed time t € (0, T) have a joint probability 
density 

Vn,t{xi, ...,Xn) = det K(aj, Xfc)] . ^^-^ „ det [pT-^ti^j^ hk)] ■ ^ , 
where Zn,t is the normalization constant such that 

/ Vn^t{xi,...,Xn)dxi---dXn = l- 

J(0,oo)" 
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This is a biorthogonal ensemble (|l.ip with functions 

fj{x) =Pt{aj,x), gj{x) =p^_^{x,bj). 

We are going to take the confluent limit Uj ^ a > 0, and bj 0. Then the biorthogonal 
ensemble structure is preserved. In our first result we identify the functions fj and gj for this 
situation. 

Proposition 2.1. In the confluent limit aj ^ a > 0, bj 0, j = 1, . . . ,n, the positions of 
the non-intersecting squared Bessel paths at time t G (0, T) are a biorthogonal ensemble with 
functions 

f2j-iix)=x^-^pf{a,x), j = l,...,ni:=\n/2], (2.4) 

f2j{x) = x^~^p"'^^{a,x), i = 1, . . . ,n2 := n - ni, (2.5) 

g,{x) = x^-'e-^ , j = l,...,n. (2.6) 

Proof. In the confluent limit Oj a, the linear space spanned by the functions y ^ pf{aj,y), 
j = 1, . . . , n, tends to the linear space spanned by 

y '-^ ^—r iP?(a,y), j = l,...,n. (2.7) 
Using the differential relations satisfied by the transition probabilities, (see e.g. [Ij or [151 [T9]): 
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g^pnx,y) = f^{prHx,y)-P?{x,y)), 



x—p'^+\x,y) = ^^Pnx,y) - (| + a + l)pr^(x,y), 

it is easily shown inductively, that the linear span of (j2.7p is the same as the linear space 
spanned by 

y ^ y^-^p^{a,y), j = l,...,ni, 

y ^ y^~^p^+^{a,y), j = l,...,n2 

which are exactly the functions in (j2.4p . (j2.5p . 

Next, the linear space spanned by the functions x i— > p^_^{x,bj), j = 1, . . . ,n, tends in 
the confluent limit bj — > to the linear space spanned by the functions 

--9^ b->T-.(^,y)]|,=o- (2-8) 

By I^A^ and we have that 

1 ^-(x+v)/(2(T-t)) 

fc=0 

which is an entire function in y of the form 



(t y) = -(x+v)/(2(T-t)) V- Wil 

y PT-tK-^^y) (2(T-t))°+i Z^yt!r(A: + a + l)(2(r-t))2fc 



fc=0 



where each Pk{x) is a polynomial in x of exact degree k. Thus the linear space spanned by the 
functions (j2.8p is equal to the linear space spanned by the functions (j2.6p . which completes 
the proof of the proposition. □ 
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Remark 2.2. In the next subsection we will see how Proposition 12.11 allows us to identify 
the ensemble of non-intersecting squared Bessel paths at any time t £ (0, T) as a multiple 
orthogonal polynomial ensemble. For the transition probability density of the (non-squared) 
Bessel process the calculations as in the proof of Proposition 12.11 would not work and in 
fact the positions of non-intersecting Bessel paths are not a multiple orthogonal polynomial 
ensemble. This is the reason why we concentrate on squared Bessel paths. 

Of course, by taking square roots we can transplant results on non-intersecting squared 
Bessel paths to non-intersecting Bessel paths, see Remark 12. 101 below. 

2.2 Multiple orthogonal polynomial ensemble 

According to Proposition 12 . 1 1 the biorthogonal ensemble in the confluent case is an example of 
what we call a multiple orthogonal polynomial ensemble. A multiple orthogonal polynomial 
ensemble in general may involve an arbitrary number of weights and an arbitrary multi- 
index, but we will discuss here the case of weight functions wqjWi, W2 and a multi-index 
(ni,n2) where ni + n2 = n and rii = \n/2\. We take functions 

/2j-i(x) = x^~'^wi{x), hjix) = x^~'^W2{x) 

and 

gj{x) = x^~'^wq{x), j = l,...,n, 

and we use these functions for a biorthogonal ensemble (jl.ip . Note that in the squared Bessel 
case, we have by Proposition 12.11 and ()2.ip that (where we drop irrelevant constants) 

«}i(x) =x'^/2g-f j^V|^^ (2.9) 

u^2(x)=x("+i)/2e-|./«+i(^) (2.10) 

wo{x) = e"^(^ (2.11) 

The biorthogonalization process leads to bases j = 1, . . . ,n, and we may take them 

so that 

(t)j{x) = Aj^i^i{x)wi{x) + Aj^i^2{x)w2{x), ipj{x) = Bj^i{x)woix), 

where and are polynomials of degrees \{j — l)/2] and [(j — l)/2j, respectively, 

and Bj^i is a monic polynomial of degree j — 1. The biorthogonality property is 

j {Aj^iwi{x) + Aj^2W2{x)) Bk{x)dx = 6j^k, j,k = 0, . . . ,n - 1, (2.12) 
where we have put 

wi{x) = wo{x)wi{x), W2{x) = 'Wo{x)id2{x)- (2.13) 

The polynomials Aj^i and Aj^2 satisfying ()2.12p are called multiple orthogonal polynomials 
of type I and the polynomials i?^ are called multiple orthogonal polynomials of type II. The 
correlation kernel 

n n—1 

Kn{x,y) = ^(l)j{x)ijj{y) = ^{Aj^iwi{x) + Aj^2W2{x)) Bj{y)wo{y) 

j=l j=0 
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is called a multiple orthogonal polynomial kernel. We will use the equivalent form (it is 
equivalent since it gives rise to the same m-point correlation functions) 



wo{y) 

which has a characterization through a RH problem, [BJ [20] 



Kn{x,y) = -^^Kn{x,y) = y^{A,^iWi{x) + A,,2W2ix)) B,{y) (2.14) 
wo[y) ^ 



Knix,y) = ^^■^l_y^ (0 wiiy) W2{y)) Y^\y)Y+{x) | 

where y is a solution of the following 3x3 matrix valued RH problem. 

1. y is analytic in C \ M. 

2. On the real axis, Y possesses continuous boundary values y+ (from the upper half plane) 
and y_ (from the lower half plane), and 

(1 'Wi{x) 'W2{xy 
1 I , x€R. (2.15) 

1 

3. Y(z) has the following behavior at infinity: 



Y{z) = { I + { - ] ] I |, z^oo, zeC\R. (2.16) 




If the weight functions are not defined on the whole real line (as it will be for the non- 
intersecting squared Bessel paths: the case of interest in this paper), we have to supplement 
the RH problem with appropriate conditions at the endpoints. The RH problem is an exten- 
sion of the RH problem for orthogonal polynomials of Fokas, Its, and Kitaev [28] to multiple 
orthogonal polynomials due to Van Assche et al. [52j . 

In this paper we have by ([2J]), ([2T0]) . ([2lT]) . and ([233]) 



wi{x) = x"/^ exp 



Tx \ r f V'^ 



2t(T -t)J"\t, , , 

V ;/ \ / ^2.17) 



..W^.<-./^e.p(-^).„«(^). 

The weights are defined on [0, cxd) so that the jump condition (j2.15p only holds for x G 1R+, 
and the RH problem ()2.15|] . (I2.16|] is supplemented with the following endpoint condition. 



4. Y{z) has the following behavior near the origin, as 2^0, zGC\ 



'1 h{z) l\ ( if -l<a<0, 

Y{z)=0\l h{z) 1 , with /i(z) = <^ log|z|, if a = 0, (2.18) 



^1 h{z) 1/ I 1. if 0<a. 

The O condition in (I2.18P is to be taken entrywise. 
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2.3 Multiple orthogonal polynomials for modified Bessel weights 

We are fortunate that the multiple orthogonal polynomials associated with the weights (j2.17p 
were studied before by Coussement and Van Assche |181ll9j. They showed that all polynomials 
Aj^i, Aj^2 and exist so that the above RH problem has a unique solution and 

dety(z) = l, forzGC\M+. 

In addition satisfies interesting recurrence and differential relations which they were able 
to identify explicitly. 

The type II multiple orthogonal polynomials satisfy a four term recurrence relation 

xBk{x) = Bk+iix) + bkBk{x) + CkBk-iix) + dkBk-2{x) 

with recurrence coefficients that are obtained from [191 Theorem 9] after appropriate rescaling 
and identification of parameters 

bk = ^ ' + \ \ 2k + a + 1), 
Cfc = k + k{k + a), 

In addition y = Bn is a solution of the third order differential equation |19l Theorem 11] 



( {n-a-2)T a \ , nT"^ 

V4t2(T-t)2^ ^ 2t{T - t) " 4t2 J ^ " 4t2(r-i)2^ " ^ ' 



2.4 Time scaling and large n limit 

We want to analyze the kernel Kn from (|2.14p in the large n limit. To obtain interesting 
results, we make the time variable depend on the number n of paths. Hence, we rescale the 
time in an appropriate way, namely we replace the variables t and T 

t 1 
1 1— > — , T — , 
2n' 2n' 

so that < t < 1. Thus, the system of weights ()2.17p now becomes n-dependent 

„,2 / nx \ /2ny/ax\ 
wi{x) = wi^n{x) = x ' exp {-j^YZT^ J ° ( 1 ) 

W2ix) = W2,nix) = CXp (^-J^Y~~tj) '^""'"^ ( t 

Alternatively, we could have performed space scaling, putting T = 1 and replacing the position 
variable x with 2nx and the starting position a with 2na. 



. (2-20) 

2n\/ax 
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After the change of time parameters t t/(2n), T l/(2n) the differential equation 
(|2.19p tm'ns into (with x replaced by z) 



zy"'{z)+(^i2 + a)-j^^ y"iz) 



( r?z n(n — a — 2) an^\ , , , , , 

Expressing (j2.2ip in terms of the scaled logarithmic derivative C = u' /{ny) and keeping only 
the dominant terms with respect to n as n ^ oo, we arrive at the algebraic equation for 

C = C(^), 

2z ^2 ( ^ 1 o.\ . 1 



which will play a central role in what follows. By solving for z, it may be written as 

\-kC 



C(i-t(i-t)C)2' 



h = (\-t){t- a(\-t)). (2.23) 



Proposition 2.3. For every t G (0, 1) the three-sheeted Riemann surface associated with 
(I2.23|) has four branch points at 0, oo, p and q with p < q. There is a critical time 

f = « e(o,i) 

a + 1 

such that 

Case 1: for t < t* we have < p < q, 
Case 2: for t > t* we have p < < q, 
Case 3: for t = t* we have p = < q. 

Note that the three cases correspond to k < 0, k > 0, and k = 0, respectively, where k is 
the constant in (j2.23p . The proof of Proposition 12.31 follows from the discussion in Sectional 

In this paper we are going to analyze Case 1 and Case 2. In order to handle the two cases 
simultaneously, we shall denote the real branch points by p- < p+ < q, where 

P- = min(0,p), p+ = max(0,p). 

Functions defined on the Riemann surface associated with (j2.23p will play a major role 
in the steepest descent analysis. There is an alternative approach based on an equilibrium 
problem for logarithmic potentials and so-called (^-functions. We briefly outline this approach 
in the appendix of this paper. 

2.5 Statement of results 

We state our results for the kernel (12.14p . 

Kn{x,y) = ^^■^l_y^ (0 wi{y) W2{y)) Y-\y)Y+{x) | | (2.24) 
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where Y is the solution of the RH problem (j2.15p , (j2.16p , ()2.18p with weights wi and W2 as 
in ()2.20p . Note that Kn depends on a > and t G (0, 1). In the following a will be fixed. To 
indicate the dependence on t we occasionally write 

Kn{x,y) = Kn{x,y;t). 

To emphasize the dependence of the branch points on t we may write p{t), q{t), P-{t), and 

Theorem 2.4. Under the resettling described above, the following hold. 

For every t G (0, 1), the limiting mean density of the positions of the paths at time t 

p{x) = p{x;t) = lim —Kn{x,x;t) 

n—^oa n 

exists, and is supported on the interval [p^{t),q{t)] C [0,oo). The density p satisfies 

p{x) = - \lmC{x)\, p+{t) <x <q{t), (2.25) 
vr 

where = (^(x) is a non-real solution of the equation (I2.23p . 

From Theorem 12. 41 it follows that as n — > oo, the non-intersecting squared Bessel processes 
fill out a simply connected region in the tx-plane given by 

< t < 1, p+{t) < x < q{t). 

This region can be seen in Figure [H 

From the definition of and q{t) as branch points of the Riemann surface for (|2.22p 

it may be shown that x = p+{t), x = q(t) are solutions of the algebraic equation 

4ax^ + _ 20at{l - t) - 8a^{l - tf) - 4x(l - t){t - a(l - t)f = 0. (2.26) 

The locus of this algebraic curve inO<t<l,x>0 gives us the boundary curve. Observe 
that it depends only on a, and is independent from the parameter a. 

There are some peculiar features of the boundary curve, which may be checked by direct 
calculation that we leave to the reader. 

Corollary 2.5. For every a > we have the following. 

(a) The lower boundary curve x = p+{t) is positive for t < t* = a/{a + 1) and it is zero for 
t > t*. At t = t* it has continuous first and second order derivatives. 

(b) The upper boundary curve x = q{t) has a slope 

q\l) = -4 

at t = 1 which is independent of the value of a. 

(c) The upper boundary curve x = q{t) is concave if a < 1. It is not concave on the full 
interval [0, 1] if a > 1. 

(d) The maximum of the upper boundary curve x = q{t) is a + 1. 
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Figure 1: Numerical simulation of 50 rescaled non-intersecting BESQ^ with a = 1 (top) and 
a = 5 (bottom). Bold line is the boundary of the domain described in Theorem I2.4[ 

By continuity the results of Theorem 12.41 and Corollary 12.51 continue to hold for a = 0, 
which is the case of non-intersecting squared Bessel bridges |35j . 

Remark 2.6. The numerical experiments leading to Figure[T]have been carried out exploiting 
the connection of the non-intersecting squared Bessel paths with the matrix- valued Laguerre 
process, as described in [361 [39]. Indeed, let a £ NU{0} and bjk, bjk, 1 < J < n + a, 1 < k < n, 
be independent one-dimensional standard Brownian motions. Consider the {n-\-a) xn matrix- 
valued process M{t) = {rrijk) with entries mjk{t) = hjk{t) + ihjkit) and define the n x n 
symmetric positive definite matrix-valued process, called the Laguerre process, by 

E{t) = M{tyM{t), te[0,+oo), 

where M(t)* denotes the conjugate transpose of M(t). Then the process of eigenvalues of 
S(t) and the noncolliding n-particle system of BESQ*^, with d = 2(a + 1), are equivalent in 
distribution. 
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Finally, in the non-critical case t ^ t* we find the usual scaling limits from random matrix 
theory, namely the sine, Airy, and Bessel kernels. 

Theorem 2.7. Let t ^ t* . Then for x* G {p+{t),q(t)), we have 

^ ( * X ^ y \ siuTrfx — y) 



n-^oo np{x*) \ np{x*)' np{x*) J tt{x — y) 
uniformly for x and y in compact subsets o/M. 
Theorem 2.8. Let t ^ t* . Then for some constant c > 0, 

1 / , , X y \ Ai(x) Ai'fw) - Ai'(x) Ai(w) 



If t < t* , then for some constant c > 0, 

^ f ^ m y ^ Ai(x)Ai'(y)- Ai'(x)Ai(y) 



n^oo cn^l'^ V cn2/3 ' cn2/3 / x — y 

Theorem 2.9. Zei t > t* . Then for some constant c > 0, and x,y > 0, 



oo cn^ ^\cn^^ cn^J \x) 2{x — y) 

In the bulk we find the sine kernel, at the soft edges we find the Airy kernel, and at the 
hard edge we find the Bessel kernel of order a. Note that the factor (y/x)"^^ in the Bessel 
kernel is not important since it will not influence the determinantal correlation functions. 
This observation also explains why totally symmetric results are obtained if we reverse the 
process and study n non-intersecting BESQ'^ paths starting at the origin and ending at a 
positive value a. Indeed, (y/x)"^"^ is the only factor in the transition probabilities (12. ip that 
is not symmetric in its variables. 

At t = t* there is a transition from the Airy kernel to the Bessel kernel. This is when the 
non-intersecting squared Bessel paths first hit the hard edge. The soft-to-hard edge transition 
is different from previous ones considered in [101 [T5] . We will treat this transition in a separate 
publication. 

Observe also that neither the boundary of the domain filled by the scaled paths, nor the 
behavior in the bulk or at the soft edge depends on the parameter a related to the dimension 
d of the BESQ'^. This dependency appears only in the interaction with the hard edge at 
a; = 0. A possible interpretation may be that a is a measure for the interaction with the hard 
edge. It does not influence the global behavior as n — > oo, but only the local behavior near 0. 

Remark 2.10. By taking square roots we can transplant Theorems 1 2 . 41 and 12 . 71(2^ to the case 
of non- intersecting Bessel paths. The correlation kernel for the positions of non-intersecting 
Bessel paths, starting at y/a and ending at is 

2ViyK„(x2,y2) 

where is the kernel (j2.24p as before. It is then easy to show from Theorems 12.71 and 12.81 
that the scaling limits are again the sine kernel in the bulk and the Airy kernel at the soft 
edges. At the hard edge however. Theorem 12.91 gives the scaling limit 

'yy y/xy Ja{x)yJ'^{y) - xJ'^{x)Ja{y) 
\xJ X + y X — y 
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The proofs of Theorems 12.41 12.71 12.81 and 12.91 are given in Section [TOl They follow from 
the steepest descent analysis of the RH problem for Y. The steepest descent analysis itself 
takes most of the paper, see Sections [HHSl 

Since we will be dealing extensively with 3x3 matrices we find it useful to use the notation 
Eij to denote the 3x3 elementary matrix whose entries are all 0, except for the (i, j)-th entry, 
which is 1. Thus 

(^ij)k,l = ^i,k^j,l (2.27) 

for i, j, k, I £ {1, 2, 3}. The following properties can be easily checked and will be used without 
comment. 



Lemma 2.11. (a) For i,j,k, I £ {1,2,3}, 

EijEki = 



Eii, ifj = k, 
O, otherwise. 



(b) If c £ C and i,j G {1, 2, 3}, i ^ j, then I + cEij is invertible, and 

{I + cEij)-^ = I - cEij . 

3 First transformation of the RH problem 

We apply the Deift-Zhou method of steepest descent to the RH problem ([2T3]) . (1236]) . (|2T8]) 
with weights wi and W2 given by p.20p and with indices ni and 112 as follows: 

{n/2, if n is even, I n/2, if n is even, 

n2 = < ' 3.1) 
(n + l)/2, if n is odd, [(n-l)/2, if n is odd. 

The steepest descent analysis has certain new features that have not appeared in the literature 
before. 

A possible approach was suggested by Van Assche et al. in [52], since the system of weights 
(|2.20p is a Nikishin system [36]. This means (in this case) that 



Wl{x) y_oo X-U 

where (T„ is a discrete measure on the negative real line, see [191 Theorem 1], with masses at 
the point 

- (t j„,fc/(2ny^))2, A; = l,2,..., (3.3) 

where ja,ki k = 1,2, . . ., are the positive zeros of the Bessel function J^. The approach of [52] 
would involve a preliminary transformation 

X = y(^I-^E2^ (3.4) 

which would result in a jump condition 

X+ = X_{I + wiEi2) (3.5) 



13 



on M_|_. Since ^ has poles on the negative real line, the third column of X has poles on the 
negative real line, which could be described by residue conditions as in [5]. We might then 
continue as in [1] by turning the residue conditions into jump conditions. However we will 
not follow this approach and we will not use the transformation (j3.4p . 

Instead, our first transformation is based on the special properties of the modified Bessel 
functions. We introduce the two functions 



yi 



(z) = z("+i)/2/,+i(2Vi), y2{z) = z(-+i)/2k,+i(2^), (3.6) 



where Ka+i denotes the modified Bessel function of the second kind, see [U Section 9.6] for 
its main properties. The functions yi and 7/2 are defined and analytic in the complex plane 
with a branch cut along the negative real axis. The jumps on M_ can be computed from the 
formulas 9.6.30 and 9.6.31 of [T]. We have 

yi+(x)=e2--yi_(x), x < 0, ^^^^ 

y2+{x) = y2-{x) - i7re*"''yi_(x), x < 0. 

From the expressions for the derivatives of the modified Bessel functions, see [H formulas 
9.6.26], we deduce that 

y[{z) = z^/^U2^), y'^iz) = -z^'^ K^{2^z). (3.8) 

The relations (|3.6p and (13. Sp imply that the weights wi and W2 defined by (|2.20p can be 
expressed in terms of the function yi and its derivative y'l as 



wi{x) =T "exp ( - _ yi(r x), 

^ (3.9) 



W2{x)=T ° ^ exp iy-j^Y~t)) y^('^'^^y 
where we have put 

T = Tn = . 

t 

We also need the following wronskian relation, see formula 9.6.15 of [1], 

y^{z)y'2iz)-y[{z)y2{z)=-^, zGC\M_. (3.10) 

Now, we are in a position to define the first transformation of the RH problem (j2.15p - 
(I2.18|) . The aim of the first transformation is to modify the jump matrix in order to have 
only one remaining weight on R+, as in (j3.5p . which is also simpler than the weights wi and 
W2- Indeed, relations (j3.9p and (|3.10p allow to remove the modified Bessel functions from the 
jumps, replacing them by a simple power function. The price we have to pay for the simpler 
jump on M_|_ will be a new jump appearing on ]R„ and on two contours that are taken as 
in Figure m We take A+ as an unbounded contour in the second quadrant asymptotic to a 
ray argz = 6 for some 9 S (7r/2, tt) as z — > oo, and meeting the real axis at the point p- < 0. 
Its mirror image in the real axis is the contour A2 . The contours A^ are the boundary of 
a domain containing the interval (— cxD,p_) and we refer to this domain as the lens around 
(-00, p_). 



14 




We define for z£C\R, 












1 


;)( 












2y2{T'^z) -z-°'yi{T'^z) 





(3.11) 



where Ci is the constant matrix 





A 













1 


,4(a+l)2 
16t 


< 







1 




(I 














1 








1 

Lot 



\ /l 

- (27rr)^^/2 Q I , if n is even, 

/ \0 i(27rT)-i/2^ 
/l 

(27rT)^i/2 I , if n is odd. 

\0 i(27rT)^i/2 



(3.12) 



Note that, in view of the wronskian relation (j3.10p . the determinant of the fifth matrix in 
the right-hand side of (|3.11|) is equal to r^". Then it is easy to see that detX{z) = 1. The 
matrix X{z) is analytic in C \ M since the matrix Y{z) is analytic in C \ M+ and yi^r'^z) and 
y2{T'^z) are analytic in C \ M_. Now define 



X{z) = X{z) 
for z outside the lens around (— cxo,p_), and 

X{z)=X{z) {I^e^^^'z-'^E^^) 



(3.13) 



(3.14) 



z 

for z in the part of the lens bounded by and (— oo,p_). [Recall that Eij is used to denote 



the elementary matrix (j2.27p .] 

From (j3.9p . (I3.10p . the jump relations (j3.7p . and the fact that Y[z) is the solution of the 
RH problem (j2.15p - (|2.18p . one derives the jump relations (j3.15p - (j3.18p below. As 2 ^ 0, we 
note the following behavior 



yi{z) 



1 



r(a + i) 
1 . 



2/2(z) ~ ^r(a + l). 



r(a)- 



y2[z) ~ < 



ir(a). 



a > 0, 
ilog(z), a = 0, 

-ir(-a)z", a<0. 
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which is a consequence of the known behavior of the modified Bessel functions near 0, see 
formulas 9.6.7-9.6.9 of [1]. This shows that X(z) has the same kind of behavior as Y(z) at 
the origin. The behavior of X{z) near the origin is then also the same, except in case p_ = 
and a > 0, see ()3.14p . The result is that X{z) is the solution of the following RH problem: 

Proposition 3.1. The matrix-valued function X{z) defined by (jS.lip . (jS.lSp . and (|3.14p is 

the unique solution of the following RH problem. 

1. X{z) is analytic m C \ (M U ). 

2. X{z) possesses continuous boundary values on (MU A^) \ {0} denoted by and X_, 
where X^ and X- denote the limiting values of X{z) as z approaches the contour from 
the left and the right, according to the orientation on M and A^ as indicated in Figure\^ 
and 

(3.15) 

(3.16) 

(3.17) 
(3.18) 




X G 



XAx 



X G (— oo,p_ 



X+(x) =X_(x)(/+|xrE32) xG(p_,0), 



X+{z)=X_{z) (l + e 
3. X{z) has the following behavior near infinity: 

X{z) =[I + 



uniformly as z 









i; 


^("1)'V4 












a/2 


»] 








-(-1)V4, 



/I 






zG Af. 





1 



V*^ -2' -2. 







~n/2^~2ny/az/t 










^-n/2^2n^/az/t 



(3.19) 



\: 



4- X(z) has the same behavior as Y{z) at the origin, see (|2.18p . either if < or if 



z ^ P- =0 outside the lens around {—oo,p-). If p^ 
{—oo,p-), then 



and z ^ in the lens around 



X{z) 




(3.20) 
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Proof. Statements in items 1, 2 and 4 are proved by straightforward calculations. It only 
remains to check the asymptotic behavior at infinity given in item 3. This follows from the 
asymptotic expansions 



y'liz) 



J_A+U2V-. (. _ 4(« + 1)^-1 , (4(a + l)^-l)(4(a + 1)^-9) ( ^\ \ 



1 a_l 
lZ2 4e 



2./i/^l (4a2-l)(4a2-9) 
16^ 512z 



1 



,3/2 



2^^ 

as z — > oo, I argz| < vr, and 

_ V^,f +i,-2Vl A + 4(a + 1)^-1 (4(a + l)^-l)(4(« + l) 
^ 2 V 512z 

, 0r - 1 4a2_i (4a2 _ i)(4a2 -9) 



(3.21) 



+ 



,3/2 



1 



,3/2 



(3.22) 

as z ^ cxD, I argz| < 37r. These formulas are consequences of the corresponding asymptotic 
expansions of the modified Bessel functions, see formulas (9.7.1)-(9.7.4) of pQ. 
It follows from (13:211) and (l3:22|) that 



A{z) 



T 

as /A 



2y'^{T'z) z-^^y^ir'z) 











-27rir^" 



(3.23) 



TXTZ 



1 i 



+ 



1 -i} ' z^l"^ \i 

as z ^ oo, I argz| < vr, where Di and D2 are diagonal matrices 

J_/4(a + 1)2-1 \ 

16r V -f(4a2-l)j' 
1 /(4(a + l)2-l)(4(a + l)2-9) 



^0(T3/2g-2rv^(T3 



D2 



512r2 











-i{Aa' - l)(4a2 -9) 



and o"3 = ( Q Jj*^ ) is the third Pauli matrix. Thus 



Aiz) 



TTTZ 



1 0\ A^^ft^ 3/.) 

-tj z 



1 i 
i 1 



^00-3/2 -2tv^(T3 



(3.24) 



where ^''2 = q*)- Now commutes with z'^^/^ since both are diagonal matrices. We also 
have 



The result is that (f3r2i|) leads to 



Aiz) 



'TTT 



V2 



TTT 



1 

-i 

1 -i 




,1/2 



-J 




Z2 







.^^3/4 



+ 0(z-i) 



yf^3/4 



1 i 
i 1 



16r 
— i 



1 fl 



^oo-3/2g-2rV^o-3 



0(73/2 -2tv^(73 



(3.25) 
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as z ^ oo, I argz| < vr. 

Now if n is even we use (j2.16p with ni = 71-2 = n/2, see ()3.ip . along with (j3.23|) . (j3.25p in 
(|3.1ip - (j3.13p to find that (j3.19p holds as z ^ 00 in the region exterior to A^. The asymptotics 
is uniform in that region. 

If n is odd then ni = n/2 + 1/2 and n2 = n/2 — 1/2, see (|3.ip . Then we need to analyze 
z~'^^/'^A{z) with A given by (j3.23p . A computation similar to the one that led to (j3.25p gives 
us 

= V2^ (^4X1 °.) {I + Oiz'^)) ^"''^^ ( • 1) z"'^«/2e-2-v^-3 

and (j3.19p follows as well, taking into account the different formula (j3.12p for the case n is 
odd. 

The asymptotic formulas p.2ip are not valid uniformly up to the negative real axis. The 
special combination y\ — ^e"'^*y2 however, does have the asymptotics ()3.2ip uniformly for 
7r/2 < argz < vr and y\ + ^e~°'^*?/2 has the asymptotics (j3.2ip uniformly for — vr < argz < 
— 7r/2. This can be seen from the formulas that connect the various Bessel functions (combine 
formulas 9.1.3-4, 9.1.35, 9.6.3-4 of [I]) 

- -e-^*y2(^) = z("+^)/2^?|i(2Vie--^/2), 

y'^{z) - -e^^%{z) = z"/2/7(i)(2^e— 

y,iz) + -e-"-*y2(^) = z^'^+'^Z^ H;^l,{2^e-'/^), 
m 

y[iz) + -e— = z^/^iJ^^) (2^e-/2), 
vrz 

where Ha^ and Ha^ are the Hankel functions, and the asymptotic expansions (see [H 9.2.7- 
10]) of the Hankel functions in the upper and lower half-planes, respectively. Hence 

yi(z)T-e^""'y2(^) 
m 

1 ^,g+i^2./i^ 4(a + 1)^-1 ^ (4(a + l)2-l)(4(a + l)2-9) ^ ^ 



2^ V 16\/^ 5122 \z^/^ 



y[{z)T-e^''^%{z) 



1 - i2/i/^ 4a2_i (4a2 _ i)(4a2 _ 9) 



2^ V 16\/i 512z V^^/^ 

(3.26) 

uniformly as z — > 00 in the region bounded by and the negative real axis. Using the 
asymptotics (j3.22p and (|3.26p . and the definition (j3.14p of X{z) in the regions bounded by 
A^ and the negative real axis, we obtain by the same calculations that ()3.19p holds uniformly 
as 2 ^ 00 in these regions as well. 

This completes the proof of Proposition 13.11 □ 
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4 The Riemann surface and the second transformation of the 
RH problem 

The Riemann surface TZ for the algebraic equation (j2.22p plays an important role in the next 
transformation of the RH problem. We repeat it here in the form (|2.23p 

-^"^^ k = {l-t){t-a{l-t)). (4.1) 



C(i-t(i-t)C)2' 

There are three inverse functions to (j4.ip . which we choose such that as z ^ oo, 

..x 1 1 t + Aajl-t) (l-t)(t + a(l-t)) / 1 \ 

1 I 1 t + Ml-O (l-t)(t + a(l-t)) , / 1 \ .44. 

Here, as in the rest of the paper, all fractional powers are taken as principal branch, that is, 
positive on with the branch cut along M_. The behavior of these functions for real values 
of z can be deduced from Figure [3] which shows the graph of z = z{C), C £M, and which also 
indicates the branches of the inverses C = Ck{z) for real z. 

It is straightforward to check that the discriminant of equation (j2.22p is equal (up to a 
non- vanishing factor depending only on t) to the polynomial in the left hand side of (|2.26p . Its 
three roots along with the point at infinity constitute the four branch points of the Riemann 
surface TZ. Analyzing the signs of the coefficients in ()2.26p it is easy to show that, according 
to the value of t G (0, 1) with respect to the critical value t = t* = a/{a + 1), the following 
cases arise (see Figure [3]): 

• Case 1: t £ (0, t*), i.e.. A; < 0. The Riemann surface TZ has three simple real branch 
points < p < q, plus a simple branch point at infinity. This is the left-most graph in 
Figure [3l 

• Case 2: t £ (t* ,1), i.e., k > 0. The Riemann surface TZ has three simple branch points 
p < < q, plus a simple branch point at infinity. This is the right-most graph in Figure 

13 

• Case 3: t = t* , i.e., k = 0. This is the critical case where the Riemann surface TZ has 
two real branch points, and q > 0, being degenerate (of order 2), and q being simple. 
The point at infinity is still a simple branch point of TZ. 

These assertions coincide with the statement of Proposition 12.31 The rest of the assertions 
of Corollary 12.51 is a consequence of straightforward although tedious computations based on 
equation (j2.26p . 

In this paper, we shall analyze Case 1 and Case 2. The sheet structure of TZ is shown in 
Figure m As before we use p- = min(p, 0) and p+ = max(p, 0). The sheets TZi and TZ2 are 
glued together along the cut Ai = q] and the sheets TZ3 and TZ2 are glued together along 
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Figure 3: Plots of z = ^(jzi(i^S)cp" ' C G IR, for Case 1 (A; < 0; left) and Case 2 {k > 0; right). 

the cut A2 = (— oo,p_]. The functions Cii C2, Cs defined and analytic on the sheets TZi, 
7^2) and TZ^ respectively, and we have the jump relations: 

Ci±{x) = (2^^), xgAi, ^^^^ 

C2±{x) = CsTi^)^ X G A2. 

On A2, the function is real and C2 and Cs are complex conjugate, while on Ai, the function 
C3 is real and Ci and (2 are complex conjugate, so that 



Ci±(2;) = C2±(a;), xgAi, 

(4.6) 

C2±(a;) = C3±(a;), x G A2. 

Near the origin, one may check from (j4.ip that, as 2: ^ 0, 

Ci(z) = i + O (z) , C2iz) = -^ + k2 + , C3(^) = ^ + k2 + (zV2) , (4.7) 

in Case 1 (p_ = 0), while 

C^(^) = -J^ + k2 + 0{z'/'), C2{z) = ^^ + k2 + 0{z^/^), Uz) = l + 0{z), 

(4.8) 

in Case 2 (p_|_ = 0), where we have set k = {1 — t){t — a(l — t)) as before, and 



t(l-t) i(l-i) 2A; 

Next, we introduce the integrals of the ^-functions. 



Ai(^) = / Ciis)ds, (4.9) 

Jq 

X2iz) = rC2{s)ds, (4.10) 

Jq 

Hz)= r C3{s)ds + X2_{P~)- (4.11) 
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Figure 4: The Riemann surface TZ, p- = and = p in Case 1, p~ = p and p+ = in Case 
2. 



The functions Ai and A2 are defined and analytic in C\ (— cxo, q], and the function A3 is defined 
and analytic in C\(— oo,p_]. By (I4.7p - (l4.8p . these functions are bounded in the neighborhood 
of each branch point p„, p^, q. By (j4.6p . 



Ai±(x) = A2±(x), xgAi. (4.12) 
From (I4.2p ~^ ()4.4p . it follows that, as z ^ 00, 



M.) = log. + ..- "-"" + °"-*» +o(l|. (4.13) 

^ ^ ^ z 1 

'*' = j(r-^) — 1 — 2'°!^'+'^ 

^t + M^^ (l-0« + a(l-<)) ^^ (4.14) 

t + 4„(l-t)^(l-0(, + a(l-<))_^^ (4.15) 



4x/^ 2z \ z3/2 



where ij, j = 1,2,3, are certain integration constants. 
We will need the following relation between £2 and £3. 

Lemma 4.1. We have £3 = £2 + tti. 
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Proof. By the definition of A2 and A3 we have for z = —R on the lower side of the cut A2, 

(A2 - A3)(Z) = r C2{s)ds - r C3is)ds = r (C2- - C3^){s)ds 

J{p-)- Jp- Jp- 

R r 

{C^_, - C^_){s)ds = / Csis)ds 

where ^r^^ is a contour from —R to p_ — e on the lower side of A2 continued with the circle 
around p- of radius e, and then back from p_ — e to — i? on the upper side of A2. Here e > 
is taken sufficiently small. Then we write 

+ / — ■^^^-zTTT^ds — [ — -ds. 

Since the integrand of the first term on the right-hand side is analytic in C \ (— oo,p_] and 
behaves as 0{s~^^'^) as s — > 00 (due to (14. 4|) ). it follows that the first term tends to as 
— > 00. For the second term we have 

^ ,ds + ^^0 



t{s - p_y/^ t 

as i? — > 00, e — > 0, and the third term is just simply — vri. Thus 



7z 



(A2 - Xs){z) = - vri + 0(1) (4.16) 

which gives the lemma in view of (|4.14p and (j4.15p . □ 
From (14. 5D. the definitions of the A-functions, and the relations 



j) Ci{s)ds = 27ri, j) (2is)ds = —27ri, 



(which follow by a residue calculation from the expansion (|4.2p of at infinity), where the 
integrals are taken on a positively oriented closed contour around Ai, we check that the 
following jump relations hold true, 

Ai+(a;) = Ai_(a;) + 27ri, x G (— cx),p+], 

Ai±(2;) = A2:f(2;), x £ Ai = [p+,q], ^^^^^ 
X2+ix) = X2~.ix) - 27ri, X G [p^,p+], 

X2+ix) = \3^{x) - 2TTi, X2-ix) = \3+{x), X G A2 = (-00, p_] . 

A straightforward consequence of these relations is the following statement: 

Lemma 4.2. Functions e^'^^^\ e^^ii^) dfid £^3(2) ^j-g analytic and single-valued outside of Ai, 
Ai U A2, and A2 respectively. 

In the sequel, we will need to compare the real parts of the A-functions on M and in 
neighborhoods of A2 and Ai. This is the aim of the next lemma. 
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Lemma 4.3. (a) The following inequalities hold true, 

ReAi < ReA2 on (0,p+) in Case 1, 
Re A3 < Re A2 on (p_ , 0) in Case 2, 
ReAi<ReA2 on (g, +00). 

(b) The open interval has a neighborhood Ui in the complex plane such that 

Re A2(z) < Re Ai(z), zeUi\{p+,q). 

(c) The open interval (— oo,p_) has a neighborhood U2 in the complex plane such that 

ReA2(-z) < ReA3(2:), 2; G C/2 \ (-oo,p_). 
The neigbhorhood U2 is unbounded and contains a full neighborhood of infinity. 

Proof. It is easy to check (see also the left-most picture in Figure [3]) that, in Case 1, C2('5) < 
Ci(s) < Csis) for s G (0,p+). Hence, from the definitions of the functions Ai and A2 and (j4.12p . 
we may conclude that ReAi < ReA2 on (0,p+). In Case 2, we have Ci{s) < Csi^) < C2{s) for 
s G (p_,0) (see right-most picture in Figure [3|). Moreover, since 



(A3-A2)(z) 



ip-h 



{C3-C2){s)ds + \2^{p^)-X2+{p^), 



we get, along with the third jump relation in (j4.17p . that Re A3 < ReA2 on (p_,0). Finally, 
on (g, 00), Ci(s) < C2{s) < Csi^) so that the third inequality in assertion (a) is indeed satisfied. 

On the -|- side of Ai, (A2 — Ai)+ is purely imaginary. Its derivative (C2 — Ci)+{^) is 
purely imaginary as well. By inspection of the Riemann surface TZ, it can be shown that this 
imaginary part is actually positive. Hence by the Cauchy-Riemann equations the real part of 
(A2 — Ai)(z) decreases as z moves into the upper half-plane, so that ReX2{z) < ReAi(2;) for 
z near Ai in the upper half-plane. Similarly, IieX2{z) < ReAi(z) for z near Ai in the lower 
half-plane, which shows assertion (b). 

The proof of assertion (c) is similar. In order to see that U2 contains a full neighborhood 
of infinity, it is sufficient to use ()4.16p . where ^/a/t > 0. □ 

A consequence of Lemma 14.31 is that we may (and do) assume that contours A^, defined 
in Section [3] (and depicted in Figure [2]) meet the real line at the branch point p^ , and lie in 
the neighborhood U2 of A2 where Re (A2 — A3) < 0. 

Using functions Xj, we can now define the second transformation of the RH problem. 



Uiz) = C2X{z) 



-n\i (z) 





V 





-n{\2{z)- 





t(l-t) 






-n{\s,(z)- 



t{T=T)) 



(4.18) 



where 



Co 




if n is even, 



if n is odd, 



(4.19) 
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and L is the constant diagonal matrix 

^e^i 



L = I I . (4.20) 

e^^ 



By Lemma 14.21 the matrix- valued function U is analytic in 

Making use of the jump relations (jS.lSp - fS.lSp for X and the definition (j4.18p one easily 
gets that the following jump relations for U: 

gn{Ai_(x)-Ai+(x)) ^a^n{\i-(x)-\2+{x)) q 
U+{x) = U^{x)\ gn(A2-(x)-A2+(x)) Q 

g"(A3-(a;)-A3+(x)) 



for X S 



.gn(Ai_{a;)-Ai+(a;)) q q 

U+{x) = U^{x)\ ^n{\2-{x)-\2+{x)) 

|^|Qgn(A3_(x)-A2+(a;)) gn(A3_ (x)-A3+(x)) ^ 

for X S (p-,0) (in Case 2), 

^gn(Ai_(x)-Ai+(x)) Q Q 

U+{x) = U-{x)\ _|^|-agn(A3-(a;)-A2+(x)) 

|^|agn(A3_(x)-A2+(2:)) q 

for X G (— cxD,p_), and 

/I 

[/+(z) = ;7_(z) 1 e=t""z-°e"(^2-A3)(2;) | ^ z e 

\0 1 

Using the jump relations (j4.17p . one checks easily that the jump properties for U simplify 
to the ones stated in the following proposition with the just introduced notation. 

Proposition 4.4. The matrix-valued function U{z) defined by (|4.18p is the unique solution 
of the following RH problem. 

1. U{z) is analytic in C \ (M U Af ). 

2. U{z) possesses continuous boundary values at (IR U A^) \ {p_,0}, and 

/I 

U+{x) = U-{x) io I , X G A2 = (-oo,p_), (4.21) 

\0 

f / + x"e"(^i-^2)W£; X £ (0,p.) in Case 1, 
U^-ix) = U-(x) X < (4 22) 

[ / + |x|"e"(^3-^2)(x)^g2, X £ {p^,0) in Case 2, 

^gn{A2-Ai)+{a;) q\ 

U+{x) = U-{x) ( gn(A2-Ai)_.(x) q\ xeAi = {p+,q). (4.23) 

1/ 
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U+{x) = C/„(x) (/ + x"e"(^i-^2)(^)^i2) , X e {q,oo), 
U+{z) = U^{z) (/ + e±""*z-"e"(^2-^^)(^)S23) , ^ G Af . 



3. As z ^ oo we have 



U{z) 



i + o[- 

z 



'10 

z^l^ 
.0 z-^l^. 



(I 



J2 V2 




(4.24) 
(4.25) 







(4.26) 



4- U{z) is bounded at if P~ < 0, and has the same behavior as X(z) at the origin, see 
(|2TgD and IK20\i . 



Proof. Jumps (|4.2ip - (j4.25p are result of straightforward calculations and Lemma l4.2i 

For the proof of the asymptotic condition in item 3 we note that property (j3.19p of X and 
the asymptotic behaviors (|4.13p - ()4.15p of the Aj-functions yield 



X{z) 





V 



i + o{ - 

z 



X L- 



e 
'I 



,0 






-{-l)"/4 



/l 

71^ 71, 



'l + 0{\) 

1-7^72 +f + 0(7^) 



'10 

,0 














l + Tlk + f + 0(ti7: 



(4.27) 



as z — > 00, where 



t + 4a(l-t) c- 
ci = Tl= ' C2 



\ — nk 



4^ ' 2 

If n is even, then by Lemma l4.ll we have that L^" commutes with all matrices before it. 
The last matrix in the right-hand side of (j4.27p can be moved to the left as in the proof of 
Proposition 13.11 The result is that ()4.27p is equal to 



'10 
1 ici 
.0 1 



1 ici \ { I + 0{ - 









1/4 








^-1/4 



/I 



V2 y/2 

\' 71^ 71 




1 





as 2 ^ 00. Then (06]) follows by the definition KTE^ - Uf^^ of U. 
If n is odd, then by Lemma 14.11 we have that 



L""diag(l,l,-1) = diag (e 



-nil —nl2 p—nl-i 
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commutes with all factors before it in (j4.27p . The result now is that (j4.27p is equal to 



/I o\ . .... /I \ A \ /I \ 

-A(l + Oll)){o zV4 ^ , 

\o -i cj ^ ^^^^\o z~^/y \o j=i j=j\o z-'^/y 

as z ^ oo, and again ()4.26p follows by the definition ()4.18p - (|4.20p of U. 

The behavior of U at the origin given in item 4 follows from the corresponding behavior 
of X, and the fact that the Xj functions all remain bounded near the origin. □ 

It follows from Lemma 14.31 that the jump matrices in (I4.22p . ()4.24p . and (j4.25p tend to 
the identity matrix / as n — > oo at an exponential rate. Moreover, by (|4.17p . (A2 — Ai)+ = 
— (A2 — Ai)_ is purely imaginary on Ai, so that the first two diagonal elements of the jump 
matrices in (|4.23p are oscillatory. In the third transformation we open a lens around Ai and 
we turn the oscillatory entries into exponentially small entries. 

5 Third transformation of the RH problem 

Here, the goal is to transform the oscillatory diagonal terms in the jump matrices on Ai into 
exponentially small off-diagonal terms. This we do by opening a lens around Ai, see Figure 
[5l We assume that the lens is contained in f/i, see Lemma 14.31 




Figure 5: Deformation of contours around Ai. 
We use the following factorizations of the 2x2 non-trivial block of the jump matrix in 

I gn(A2-Ai)_(x) ) 

_ / 1 OA / / 1 OA 

We set 

for z in the domain bounded by A]*^ and Ai, and we let 

T{z) = U{z) (5.2) 
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for z outside of the lens. 

Straightforward calculations show that T{z) is a solution of the following Riemann-Hilbert 
problem which is stated in the next proposition. 

Proposition 5.1. The matrix-valued function T{z) is the unique solution of the following 
RH problem. 

1. T is analytic m C \ (M U A]*" U A^). 



2. T has a jump Tj^{z) = T^{z) jxiz) on each of the oriented contours shown in FigurelM 
They are given by 







3t{x) 
3t{x) 



\x\ 




X G A2, 



ze A± 



'1 


r / + x°e'^(^i-^2)(x)£;^2, X G (0,p+) in Case 1, 

\ / + |x|"e"(^^--^2){^)^32^ J. £ (p_,o) in Case 2, 



X G Ai, 
































jj,{z) = 1 + z-"e"(^2-Ai)(^)_g2^^ ze Af, 
j^{x) = 1 + x"e"(^i-^2)(x)^^2^ ^ g 

3. As z ^ 00, we have 
T{z) = 









(1 





^1/4 







i 







1° 






1 





4- For — 1 < a < 0, T{z) behaves near the origin like: 

T{z)- 

For a = 0, T{z) behaves near the origin like: 




as z 



0. 



T{z) 




as z 



outside the lens that ends in 0, 



and 



T(z) 



O 



O 



/l logl^l log|z| 
1 log|2;| log|z| 

\1 log|z| log|z| 

^log|z| logl^l 1^ 
log|2:| log|2:| 1 

\log|2;| log|2;| 



(5.3) 



as 2; — > inside the lens around A2 in Case 1, 



as z 



inside the lens around Ai in Case 2. 
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For a > 0, T{z) behaves near the origin like: 



'I 1 l^ 

T{z) =0\l 1 1 
,1 1 l> 



and 



T{z) 



1 


1 


\z 


—a 


1 


1 


\z 


—a 


1 


1 


\z 


—a 


\z 


—a 


1 




\z 


—a 


1 




\z 


—a 


1 





as z — > outside the lens that ends in 0, 



as z — > inside the lens around A2 in Case 1, 
as z — > inside the lens around Ai m Case 2. 



5. T is hounded at p and q. 

Proof. All properties follow by straightforward calculations. 

Because of the prescribed behavior at the origin, it is not immediate that the RH problems 
for U and T are equivalent. Reasoning as in [41^ Lemma 4.1] we can still show that they are. 
Thus in particular the solution of the RH problem for T is unique. □ 



6 Model RH problem for the global parametrix 

In view of Lemma 14.31 the jump matrices in the RH problem for T all tend to the identity 
matrix exponentially fast as n — > 00, except for the jump matrices on Ai and A2. Thus we 
expect that the main contribution to the asymptotic behavior of T is described by a solution 
Na of the following model RH problem. 

1. Na is analytic in C \ (Ai U A2). 

2. Na has continuous boundary values on Ai U A2, satisfying the following jump relations: 



xGAi, (6.1) 
, 2; G A2. (6.2) 









( 





x" 





iV„+(x) 


= Na- 


(x) 












































Na+{x) 


= Na- 


(x) 










X 










1x1° 








3. As z — > 00, z G C \ A2, 



Na{z)= [l+0[- 















^1/4 











^-1/4 







-a/2 



(6.3) 
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The asymptotic condition at infinity looks a bit awkward. However it is consistent with 
the jump on A2 since it may be checked that 

satisfies 

Sh.(x)=B_(x)(|^°„ ~'o")' 

We solve the RH problem for in two steps. First we solve it for the special value a = 
and then we use this to solve it for general values of a. In both steps we will use the mapping 
function (jlTT]) 

1-kC 

^ " c(i - ccr 

with 

c = t{l-t) and k = {l-t){t-a{l-t)), (6.4) 

which gives a bijection between the Riemann surface TZ and the extended ^-plane. The 
mapping properties are summarized in Figure [6] for the two cases (Case 1 in the upper part 
and Case 2 in the lower part of the figure). 
The figure shows the domains 

nj = Cj{n,), i = 1,2,3, 

where TZj is the jth sheet of the Riemann surface, and also the location of the points 

1 1 



Cp = C2(p), Cq = C2iq), Coo=C2(oo) 



t(l-t) 



for the two cases. We observe that (2+{^i) and C2+{^2) are in the upper half plane, while 
C2-(Ai) and (^2~i^2) are in the lower half plane. 

To solve the model RH problem for a = 0, we use the polynomial -D(C) 

^(C) = (C-Cp)(C-C,)(C-Coo). (6.5) 

The square root Z)(C)^''^, which branches at these three points, is defined with a cut on 
C2-(Ai) U C2-(A2), which, as noted before, are the parts of the boundary of 7^2 that are in 
the lower half of the ^-plane. We assume that the square root is positive for large positive C- 

Proposition 6.1. A solution of the model RH problem for Nq is given by 

/Fi(Ci(z)) Fi(C2(^)) F^{Uz))\ 
No{z) = F2{Ci{z)) F2{Uz)) F2iC3{z)) (6.6) 

Vi^3(Cl(^)) Fs{C2iz)) Fs{C3iz))J 

where 

^l(0 = ^i%g?#, ^^(0=i^2^^, ^3(0 = /^3%(=^, (6.7) 

with -D(C) given by (j6.5p . Furthermore, C* 7^ Coo, cLnd Ki, K2, are explicitly computable 
non-zero constants that depend on a and t. 
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Ai 



p q 





^2 



L Co 



•C2-(Al) 




C2-(A2) 



Figure 6: Bijection (j4.ip between the Riemann surface TZ and the extended ^-plane in the 
Case 1 (top) and 2 (bottom). 



Proof. Note that each of the functions Fj, j = 1,2,3, defined in (|6.7|) satisfies 

fF,+ (C) = -Fj.iC), C G 5^2 n {ImC < 0}, 
[Fj+iO = Fj.iC), C G 57^2 n {ImC > 0}, 

because of the choice of the branch cut for D((^)^^^. From this it follows that the jth row 
[Nji, Nj2, Njs) of A'^o given in (|6.6p has the following jumps on Ai: 

(iV,i)+(z) = -{Nj2)-{z), 
(iV,2)+(z) = (iV,i)_(z), ZGAi, 
(iV,3)+(z) = (iV,3)-(^), 
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and the following jumps on A2: 

(iV,i)+(z) = (iV,i)_(z), 

(iVj2)+(z) = {Nj2)^{z), ZGA2. 

(iV,3)+(z) = -{N,2)-{z), 

These are exactly the jumps required by ()6.ip and ()6.2p when a = 0. 

It remains to verify the asymptotic condition (j6.3p with a = 0. Since the computations 
are straightforward but cumbersome, we give here the outline of the argument. Observe that 

Ci(oo) = , C2(oo) = C3(oo) = Coo ■ (6.8) 

Function Fi verifies 

Taking into account (j6.8p and (j4.2p - (j4.4p . we get that as z — > 00, 

iVn(^) = i^i(Ci(^)) = K^ + 0{l/z) , 

Nuiz) = Fi(C2(z)) = O ^ ^^3(^) = Fi(C3(z)) = O . 

With Ki = C^2/)(o)i/2 it yields 

iVn(z) = 1 + 0{l/z), Ni2{z) = ^ 



as 2 — > 00, which matches the asymptotic condition for the first row of Nq in (j6.3p . 
Analogously, 

F2(c) = o(C), C^o, 

i^2(C) =PliC- Coo)"^/' + /32(C - Coo)^/' + O ((C - Coo)=^/') , C - Coo , 

where f3i and /32 are explicitly computable in terms of K2, C* &iid the rest of the parameters 
of TZ. By (j4.2p - (j4.4p . and taking into account the second relation in (j4.6p . we have 

N2i{z)=F2{Ci{z))=0{l/z), z^oo, 

iV22(^) = F2{C2{Z)) = Z^l^ (^1 + %Z-^I'' + ^3^-' + O (^^'/')) , ^ ^ 00 , 
iV23(^)=i^2(C3(^))=i^'/' (^1- ^2^-1/2 +^3Z-l + 0(z-3/2^) , Z 



00 



where again /3j's are explicit. Imposing the condition that /3i = l/\/2 and /32 = 0, which 
determines K2 and C*; we obtain that for a certain constant 02, 

iV2i(z) = 0(l/z), 

iV22(.) = ^.^/^(l + f + 



iV23(z) = ^zl/Ml + ^ + 0(. 
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matching the asymptotic condition for the second row of A'^o in (j6.3p . 
Finally, 

m) = 7i(c - Coo)'/' + o({c- Coo)'/') , c - Coo , 



where 71 is explicitly computable in terms of and the rest of the parameters of TZ. By 
(|4.2p - (j4.4p . and taking again into account the second relation in (j4.6p . we have 

N3i{z) = Fs{Ci{z)) = 0{l/z), z^oo, 

Nssiz) = FsiCsiz)) = -i z'/^ (71 - 72^-'/' + O , z ^ oo , 

where again 7j's are explicit. Imposing the condition that 71 = i/\/2> which determines K^, 
we obtain that for a certain constant 03, 

N3i{z)=0{l/z), 

^32(.) = ^.-'/^ {l + ^ + Oiz-^)), 
A^33(.) = ^.-/^(l-^ + 0(.-)), 

as z ^ 00. This is precisely the asymptotic condition for the third row of Nq in (j6.3p . This 
concludes the proof. □ 

To construct the solution for general a, we use functions 

ri(C) =log(l-cC), CG^i, 

r2(C) = log(l-A:C)-logC-log(l-cC), C e ^2, (6-9) 

r3(C) =log(l-cC)+i7r, C e ^3- 



where c and k are as in (j6.4p . The branches of the logarithm are defined as follows. 

• log(l — cC) vanishes for C = 0> ^-nd has a branch cut along C2-('^2) in Case 1, and along 
C2-(A2) U [Cp, +00) in Case 2. 

• log(l — A;C) vanishes for = 0, and has a branch cut along (—00, 1/k] in case 1 (when 
A; < 0), and along [l/k, +00) in Case 2 (when k > 0). 

• logC is the principal branch of the logarithm, i.e., with a cut along (— oo,0]. 

With these definitions of the branches we have that rj is defined and analytic on TZj for 
j = 1, 2, 3. To see this for j = 2, it is important to note that l/k is in TZi in Case 1, and l/k 
is in 7^3 in Case 2. 
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Proposition 6.2. A solution of the model RH problem for general a is given by 

N^iz) = Co,No{z) { e"^2W q (g.io) 

gaGsiz) I 




where Nq is given by (|6.6p . 

Gj{z) = rj{Cj{z)), j = 1,2,3, zGTZj, (6.11) 
with ri, r2, r^ defined in (|6.9|) . and Ca is a constant matrix given explicitly in (|6.14|) below. 

Proof. From the definitions (j6.9p with the specified branches of the logarithm it follows that 
the functions r^, j = 1, 2, 3, satisfy the following boundary conditions 

r2(C)=n(C) + logz, Ce5^i, 
^^2(0 = r3(C) + logkl, Cg5^3, 

where z = is given by ()4.ip . Then by (j6.12p and (j6.1ip we obtain 

G2±{x) = \ogx + Gi^{x), xEAi, ^^^^^ 
G2±(2;) = log|x| + G3=p(x), x E A2. 

Using (j6.13p and the jump propertes (j6.ip . (j6.2p for a = 0, it is then an easy calculation to 
show that (j6.10p satisfies the jump conditions ()6.2p and (|6.ip . 
We note that by ^Mi and (l6lTD . 

gGiW = l_cCi(z), e^2(^) = z(l - cC2(^)), e^3(^) = cC3(^) - 1, 
and thus as z ^ 00 by ([42]), (1131), dilill . 

gGi(.) = 1 + 0(1/^)^ 

,G.(.) ^ ,1/2 - t) + ^ + t[l-t)[t^^ail-t)) ^ ^^^_3/,, 

y 1^ z olz 

y ZyJ Z 0\J CLZ 

To obtain ()6.3p we should then take the constant prefactor Ca in ()6.10p as 
/I /I 

c„= V^(i-t) = (V^(i-t))-" 

Vo V^{i-t)J \o (V^(i-t))-° 

(6.14) 

Then with the choice of ()6.14p , we indeed have that Na defined in (IG.lOp satisfies the conditions 
in the model RH problem for general a. □ 

Lemma 6.3. The solution of the model RH problem given in Proposition \6.S\ has the 
following behavior near the branch points 
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(a) In Case 1 we have 



O 



\z 


- Q 


-1/4 


\z 


- Q 


-1/4 


\z 


- Q 


-1/4 


\z 


- Q 


-1/4 


\z 


- q 


-1/4 


\z 


- q 


-1/4 


\z 


-P 


-1/4 


\z 


-P 


-1/4 


\z 


-P 


-1/4 


\z 


-P 


-1/4 


\z 


-P 


-1/4 


\z 


-P 


-1/4 



as z 



as z 



P, 



(6.15) 



(6.16) 



and 



-a/2 




o 



1/4 


\z 


-1/4 


1/4 


\z 


-1/4 


1/4 


\z 


-1/4 



as z 



0, 



(6.17) 



(b) In Case 2 we have 



and 



N^{z) = O 



\z 




|-l/4 


\z 


- q 


|-l/4 




\z 


- ^ 


|-l/4 


\z 


- Q 


|-l/4 




\z 


- ^ 


|-l/4 


\z 


- 1 


|-l/4 




1 


\z 


-p\~ 


1/4 


\z 


-P\~ 


1/4 


1 


\z 


-p\~ 


1/4 


\z 


-P\~ 


1/4 


1 


\z 


-p\~ 


1/4 


\z 


-P\~ 


1/4 



as z 



as z 



P, 



(6.18) 



(6.19) 



0> 

Na{z) I 

L 



z 


-1/4 


\z 


-1/4 




z 


-1/4 


\z 


-1/4 




z 


-1/4 


\z 


-1/4 





as z 



(6.20) 



Proof. Observe that for j = 1, 2, 3, Fj{C) = O {{( - Cg)"^-^^) as C ^ Cq, where F/s are defined 
in (j6.7p . Furthermore, for the mapping (j4.ip . (^^{(g) = (^^{(q) = q, and C3~^(C(?) is a regular 
point of TZ. Since functions Ci and ^2 are bounded and have a square root branch at q, by 
definition (16. 6|) we obtain ()6.15p for Nq. Since the transformation in (I6.10p does not affect 
the behavior at q, this proves the first identity of the Lemma. The rest of the conditions is 
analyzed in a similar fashion, and we omit the details. □ 

By (I6.ip - (16.2j) . we have that detA^^c^ is analytic in C \ {0,p,q}, and by Lemma 16.3^ 
det Na{z) = O {\z — zqI"^/^) as z — > where zq is any one of the branch points 0, p and q. 
Hence, detA^a is entire. From (j6.3p . lim^^^oo det A''q(z) = 1, and we conclude that 

detiVc,(z) = 1, zgC. (6.21) 

Comparing the local behavior in Proposition 15. II and Lemma [6. 3 1 we see that near the branch 
points, the matrix TN~^ is not bounded which means that A'^^ is not a good approximation 
to T. Hence we need a local analysis around these points. 



34 



7 Parametrices near the branch points p and q (soft edges) 



We are going to construct a local parametrix P around q. The local parametrix around p can 
be built in a similar way, and is not further discussed here. Consider a small fixed disk Bs 
with radius S > and center at q that does not contain any other branch point. We look for 
a 3 X 3 matrix valued function P such that 




Figure 7: Construction of a parametrix around q. 



1. P is analytic in 5^ \ (M U Af ). 

2. P has a jump P+{z) = P^(z) jxiz) on each of the oriented contours shown in Figure [3 
given by the restriction of Jt in Proposition 15.11 to these contours. Namely, 

/ x° 0\ 
jrix) = , x£{q-6,q) = AinBs, 

V 1/ 

jriz) = 1 + z-"e"(^2-^i)(^)E2i, zGAfn Bg, 

j^(x) = / + x°e"(^i-^2)(^)£;i2, xe{q,q + 6). 



3. As n — > oo, 

P{z) = Na{z){I + 0(1/?!)) uniformly for z e dBs \ {RU Af), 
where A^^ is the global parametrix built in Section [6l 

4. P is bounded as z ^ q, z £ M. \ Af. 

The solution of the RH problem 1.-4.) can be built in a standard way using the Airy 
functions; we follow the scheme proposed in [22l [231 125] and developed, for instance, in 
[3 [211132]. The function 



-(A2-Ai)(z) 



2/3 



(7.1) 
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is a biholomorphic (conformal) map of a neighborhood of q onto a neighborhood of the origin 
such that f{z) is real and positive for z > q. We may deform the contours Af near q in such 
a way that / maps Af D Bs to the rays with angles ^ and — respectively. We put 

yo{s)=Ai{s), yi{s) =ujAi{ius), y2is) = cu^ Ai{LU^ s) , uj = e^""'/^ , 
where Ai is the usual Airy function. Define the matrix ^ by 

^yo(s) -y2{s) 0^ 
^(s) = I y'ois) -y'2{s) 

1 

-yi{s) -y2{s) 
^(^) = I -y'lis) -y'2{s) 

1 

-y2{s) yi{s) 0^ 
^(^) = I -y2(s) 

1 

/yo(s) yiis) ON 

^(^) = yo(5) y'M 
V 1, 

Then (see e.g. (22l Section 7.6]), for any analytic prefactor E, we have that 

P{z) = E{z)^ (n2/3 /(^)) diag (^^-"/2et(A2-Ai)(.)^ ^„/2^-5(a,-A0{.)^ 1^ (7.2) 

satisfies the parts 1.-3. of the RH problem for P. The freedom in E can be used to satisfy also 
the matching condition (4) . The construction of E uses the asymptotics of the Airy function 
Ai(s) as s — > oo, and follows the scheme, exposed in the literature (see e.g. [40]), and we omit 
the details here. The result is the following. 

Proposition 7.1. The matrix-valued function P given in (j7.2p with E given by 
E{z) = iV,(z) 

V l/ 

0\ (n^l^f^lHz) 
i^AF n~'^/^r^/\z) I , (7.3) 

1/ V 

satisfies all conditions 1.-4- the RH problem for P . 



.a/2 












•) 


















•(- 


-i^/n 










8 Parametrix near the branch point (hard edge) 

From the local behavior of T{z) as z ^ oo, described in Proposition 15.11 it follows that the 
local parametrix P at the origin will be different from the parametrices at the other branch 
points. Fortunately, this kind of behavior has been analyzed in [IT] (for a 2 x 2 matrix valued 
RH problem) , and in [35] (for a 3 x 3 matrix valued RH problem) and we use the construction 
from these papers. There will be a new feature though in the case — 1 < a < 0. 
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8.1 Case 1. 



Let Bs be a small fixed disk with radius 5 > 0, now centered at the origin, that does not contain 
any other branch point. Consider all the jumps matrices Jt of matrix T on curves meeting 
at 0, see item 2. in Proposition 15.11 The off'-diagonal entry in jj- on (0,6) is x°'e^^^'^~^^^^^\ 
which is exponentially small since Re (Ai — A2) < — c < on [0, 6). This suggests that we may 
ignore the jump on (0, 6) in the construction of the local parametrix. Note however, that 
is not bounded as x — > in case —1 < a < 0, which explains why we need an extra argument 
for this case. 




Figure 8: Contours for the local parametrix around in the Case 1, for a > (left picture) 
and for —1 < a < (right picture). 



8.1.1 First part of the construction (which works for a > 0) 

In this part we simply disregard the jump matrix on (0, 6). Taking into account Proposition 
15. H we thus look for a 3 x 3 matrix valued function Q such that 

1. Q is analytic in Bs \ (A2 U ). 

2. Q has a jump Q+{z) = Q^{z)jT{z) on (A2 U A^) n i?^, see the left picture in Figure [HI 
The jump matrices are given by 

jt(x)=|o I , xe (-5,0) = A2nB5, 






j^{z) =1 + e^"^^z-"e"(^2-^=')(^)E23, z G Af n Bs- 

3. For — 1 < a < 0, (5(z) behaves near the origin like: 

/I \zf 1^ 

Q{z) = o[l \zf 1|, asz^O. (8.1) 
Vl Ul" 1 
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For Q = 0, Q{z) behaves near the origin hke: 

(1 logl^l logl^l 
1 log|z| log|2;| I , as z ^ 0. (8.2) 
1 log |z| logl^l 

For Q > 0, Q{z) behaves near the origin like: 

/I 1 

Q{z) = 011 1 \z\~'^ , as z — > in the lens around A2, bounded by A^, (8.3) 
\1 1 

/I 1 1\ 

Q(z) = 01 1 1, asz^O outside the lens. (8.4) 

V 1 1/ 

4. As n ^ 00, 

Q(z) = Na{z){I + 0{l/n)) uniformly for z£dBs\ (A2 U Af), (8.5) 
where Na is the parametrix built in Section [6j 
Consider 

Q{z) = Q{z) diag (1, (±l)"2-"/2ef (A2-A3){.)^ (±l)n^a/2g-f (A2-A3)(.)^ ^ fo^, ^ > 

(8.6) 

where z"/^ denotes the principal branch, as usual. By Lemma [4. 21 the diagonal factor in ()8.6p 
is analytic in Bs \ {—S, 0). It follows that the matrix valued function Q should satisfy: 

1. Q is analytic in \ (M U Af ). 

2. Q has a jump Q+{z) = Q^{z)jQ{z) on each of the oriented contours shown in Figure [HI 
left. They are given by 

/I 0\ 

j^{x)=[0 -1 , x£{-6,0) = A2nBs, 

Vo 1 0/ 

j^iz) = 1 + e±°-^^23, zeAfn Bs, 
where we have used the last identity in (j4.17p . 

3. For —1 < a < 0, Q{z) behaves near the origin like: 

/I |z|"/2 |z|"/2\ 

Q{z) = [1 |z|"/2 |z|"/2 1 as z ^ 0. (8.7) 

\i \z\"/^ \z\"/y 

For a = 0, Q{z) behaves near the origin like: 

(1 logl^l log|z| 
1 logjzl log|z| I , as z ^ 0. (8.8) 
1 log Izl loglzl 
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For Q > 0, Q{z) behaves near the origin hke: 
Q{z) = O 



and 



Q{z) = O 



/2 


\z\~ 


a/1 


/2 


\z\~ 


a/2 


/2 


\z\^ 


a/2 


a/2 


\z 





as z — > outside A 



-a/2 |2|-"/2 
-a/2 I2I-Q/2 



as Z 



2 ) 



inside A^. 



(8.9) 



(8.10) 



Although we have different expressions for Q for the cases n even and n odd, there is no 
distinction between these two cases in the conditions on Q. 

The problem for Q has a solution in terms of the modified Bessel functions of order a 
see |41l Section 6]. Namely, with the modified Bessel functions and K^, and the Hankel 
functions and H^^ (see [H Chapter 9]), we define a 2 x 2 matrix ^(C) for | argC| < 27r/3 
as 

' UK'") ^K^e") \ 

,27r<i/2j;^(2^i/2) _2CV2k^(2C1/2)J ' 
For 27r /3 < arg < vr we define it as 



*(C) 



.11 



ii7«(2(-C)^/2) 



^(C) 

And finally for — vr < argC < — 27r/3 it is defined as 

^(0 = 



ii7(^)(2(-C)V2) 



^vrC^^ /^(DV (2(_c)i/2) ^c^/2 (^(2)y (2(-C)V2) 



(8.12) 



\if(^)(2(-C)V2) 



iif«(2(-C)V2) 



f//i'))'(2(-C)V2) vrCV2 (//^V (2(-C)V2)^ 



e 2 



Q7r2cr3 



.13) 



Then we define a 3 x 3 matrix ^, given in block form by 



^(C) 







(ji^(C)ai 



0-1 



1 

1 



.14) 



[The conjugation by ai is needed to interchange the second and third rows and columns.] 
The function 



- 2 






.2 L 



can be continued analytically from Bs \ {—6,0] to the full neighborhood Bs, giving a biholo- 
morphic (conformal) homeomorphism of a neighborhood of the origin onto itself (see (j4.7p ) 
such that f{x) is real and positive for x E (0,5). Again, we may deform the contours A^ 
near in such a way that / maps A^ n Bs to the rays with angles % and — respectively. 



It follows from 



that for any analytic prefactor E, we have that 



Qiz)=E{z)^f{n^f{z)) 
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satisfies the conditions (1), (2), (3) needed for Q. So we complete the construction of Q by 
defining 

Q{z) = E{z)^{n^ f{z))dmg{l, ^^/Sg-f (Aa-As)^^ ^-«/2gt(A2-A3){^)^ ^ (8_;L5) 

where analytic in B^, is chosen to satisfy the matching condition on dB^. Using again the 
results of [H] , and taking into account that we have to interchange the second and third rows 
and columns, we define 

=iV„(z)diag(l, z"/2)diag(l,i= ~^y^ 

xdiagfl, (2T:n)-^/^ f{z)-^'\ {2T:nf/^ f{zY'^) ■ (8.16) 



Here the branch of f^/^{z) is positive for z G (0, 6). Observe that f^^^{z) = O (-z^^^) as z — > 0, 
so by (|UT7D . 

/I 1 1\ 

E{z) = 01 i\ as z ^ 0. 

It is easy to check that 

/I OA 

^+(x) = E_(x) , x€i-6,0). 

\0 -t{U/f.Yl\z)l 

Since f]f^{x) = if^^^{x) for x G (—5, 0) and E cannot have a pole at the origin, we conclude 
that E is analytic in Bg. 

Finally, the matching condition (jS.Sp in condition (4) of the RH problem for Q is satisfied 
by results of [H]. We have thus established the following. 

Proposition 8.1. The matrix-valued function Q defined by (jS.lSp . (j8.16|) . with ^ as in (j8.14p 
satisfies the conditions 1.-4- of the RH problem for Q . 

Taking into account (j6.2ip and that det^ = 1 (see |41j ) we also conclude that 

detQ(z) = l, zGBs. (8.17) 
If we would take Q as the local parametrix for T, we would define the final transformation as 

Riz) = Tiz)Qizy\ zeBs. 
Then R would be analytic in Bg \ (0, oo) with the following jump for x G (0, 6), 
R^{xy^R+{x) = Q{x)T^{xy^T+{x)Q{xy^ 

Lemma 8.2. For a > 0, the matrix Q{x)Ei2Q{x)~^ is bounded as x — > 0, x > 0. 
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Proof. If a > 0, then it follows from (j8.4p and (j8.17p that both Q{x) and Q{x) ^ are bounded 
as X — > 0, x > 0, and the lemma follows. 

For a = 0, the above argument, now based on (|8.2|) instead of (|8.4|) . does not work, since 
it would lead to a bound ©(log |x|) as x — > 0. To prove the lemma for a = 0, we look at the 
precise construction of Q. From (j8.11|) and the known behavior of /o(C) ^^iid Kq{() as C — ^ 0, 
we obtain 

^^^> ""[l log Id 
Since det ^'(C) = 1, it then follows by (l8l^ that 

^ /I 

*(C)^^ = 0{1) 0{l) I as C ^ 0. 

\0 0(log|C|) 0(log|C|), 

Using this in ()8.15p we obtain 

/I 
Q(x)-'^=io 0(1) 0(1) \ E~'^{x), asx^O,a;>0, (8.19) 

\0 C'(log|x|) C'(log|x| 

where E~^{x) is bounded near x = 0. Since 



Q{x)Ei2Q{xy^ = Q{x 



and Q{x) is bounded by (fO) and (O 1 O) Q{x)-^ is bounded by (f8?T9]) . the lemma 

follows for a = as well. □ 

From Lemma 18.21 and the fact that Re (Ai — A2) < — c < 0, for some c > 0, it follows that 
the jump matrix (j8.18p is exponentially close to the identity matrix as ?i — > 00, uniformly for 
x € (0, (5), in case a > 0. We take the parametrix P = Q in case a > 0. 

This does not work if a < 0, since then we would get that Q{x)Ei2Q{x)~^ is of order 
x° as X — > 0. Then for any fixed x > 0, the jump matrix is close to the identity matrix as 
n — > 00, but it is not valid uniformly for x G (0,6). 

8.1.2 Second part of the construction, for — 1 < a < 

Let us analyze now the case when — 1 < a < 0. Now we cannot simply ignore the jump 
matrix of T on (0, 6), so we will try to match all four jumps. Namely, we build a 3 x 3 matrix 
valued function P such that 

1. P is analytic in Bs \ (A2 U U (0, 5)) . 

2. P has a jump P-\-{z) = P-.{z) jx^z) on each of the oriented contours shown in the right 
picture of Figure [HI The jump matrices are given by 



jT{x)=\0 -|x|-° I , xe{-S,0) = A2nBs, 
Ixl'^ 
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j^(z) =1 + e±"'^'z-"e'^(^2-^^)(")^23, z£A^n Bs, 
= / + x"e"(^i-^2)(^)^i2, X G (0, 5). 

3. P{z) behaves near the origin like: 

P(z) =01 1 , as z ^ 0. (8.20) 

\1 1/ 

4. As n — > oo, 

P(z) = Ar^(z)(/ + C)(l/n)) uniformly for z G 95^ \ (A2 U Af U (0, 5)) , (8.21) 
where A'^q, is the parametrix built in Section [6j 

We use the matrix- valued function Q given by formulas (j8.15p and (j8.16p . that worked as 
a parametrix for the case q > 0. We take P in the form 

P{z) = Qiz)S{z), (8.22) 

where S is given in the four components of Bs \ (A2 U A^ U (0, 6)) as follows: 

5(.)=I + ^-l-z"e«(^-^^)(^)i5;i2, 

for z in the region bounded by (0,6) and A^, (8.23) 

2o7ri 

for z in the region bounded by (0,5) and A^, (8.24) 

S(z) = I + L_z°e"(^i-^2){^)^ !!^e"(^i-^3)(^)^i3, 

for z in the region bounded by A2 and A^, (8.25) 

5(z) = / + , . z"e"(^^~^^)(")^i9 + ^ . e"(^^"^^^("^^i3, 

for z in the region bounded by A2 and Ag • (8.26) 

This construction is actually valid for any non-integer a. 

It is a straightforward, although somewhat lengthy, calculation to show that P satisfies 
all the jump conditions from item 2. in the RH problem for P. To check the jump on 
A2 = (—5,0) one has to keep in mind that A2+ = A3_ — 2Tri on A2, see (I4.17p . and that 2;" 
is defined with a cut on (— oo,0]. The conditions 1., 3., and 4. in the RH problem for P are 
easy to verify from the above definitions and the corresponding conditions in the RH problem 
for Q. For condition 4. we also need to note that Re (Ai — Xj){z) < — c < for j = 2, 3 and 
z G Bs. 

In order to unify notation for a > and —1 < a < 0, we take as the parametrix in Bs 
in the Case 1 the matrix valued function P = QS, where 5 = / if a > 0, and S is given by 
(I823])-(I826]), if -1< a < 0. 
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8.2 Case 2. 



The construction of the local parametrix P near the origin in Case 2 follows along similar lines 
as the construction in Case 1. In Case 2 the geometry of the curves in the RH problem for 
T is shown in the right picture of Figure [H Now the jump matrix on {—S, 0) is exponentially 
close to the identity matrix if n is large, and in the first step of the construction we ignore 
the jump on (—5,0), thereby giving us the contours as in the left picture of Figure [D 




Figure 9: Contours for the local parametrix around in the Case 2, for a > (left picture) 
and for — 1 < a < (right picture). 



8.2.1 Construction for q > 

We start by constructing a solution to the following RH problem (see left picture of Figure 
i). 

1. Q is analytic in Bs \ (Ai U Af ). 

2. Q has a jump Q+{z) = Q-{z) jriz) on each of the oriented contours shown in Figure 
[9l They are given by 

jrix) 





( 






— X 


Q 









= I + z- 





3. For —1 < a < 0, Q{z) behaves near the origin like: 

/I 1\ 

Q{z) = 01 |z|" 1 , as z ^ 0. (8.27) 
\1 1/ 

For a = 0, Q{z) behaves near the origin like: 

(log|z| log|z| l\ 
log|z| log|z| 1 , as z ^ 0, (8.28) 
log|z| log|z| 1/ 
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For < a, Q{z) behaves near the origin hke: 

as z ^ in the lens around Ai, 



1 r 

Q{z) = o I 1 1 

1 1, 



/I 1 i\ 

Q{z) = 011 1 ll, asz^O outside the lens. 

V 1 1/ 

4. As n ^ oo, 

Q(^) = Na{z){I + ©(I/Ji)) uniformly for zedBs\ (Ai U Af ) 
where Na is the parametrix built in Section [6j 
With ^ built in ()8.1ip - (j8.13p we define a 3 x 3 matrix-valued function 



^.29) 
^.30) 



(8.31) 



V 

where now ^ is in the upper left block, and 



0-3 



1 
-1 



(8.32) 



Q(z) = £;(z)$(n2/(^))diag((±l)"(-z)-"/2et(A2-Ai){^)^ (±l)n(_^)"/2g-f (A2-Ai){^)^ ^ 

(8.33) 

for ±lmz > 0, where (— z)"^/^ is positive for z G {—6, 0) and is defined with a cut on (0, +oo). 
Here / is the conformal map 



i(A2-Ai)(z)-i(A2-Ai)(0) 



- 2 


- /■ 




.2 Jo 



(C2-Ci)(s)(is 



.34) 



and the analytic prefactor E is 



1 /I i 



£(z) = N^{z) diag ((-z)"/2^ (-^)-"/^ l) diag ^ . , 

X diag ((2^n)i/2/(z)i/4^ (27rn)-i/2/(^)-i/4^ 1^ ^ (8 35) 

Then we find the following analogue of Proposition 18. 1[ 



Proposition 8.3. The matrix-valued function Q defined by (j8.33p . (j8.35p . with ^ as in (j8.32p 
and f as in (I8.34p . satisfies the conditions 1.-4- of the RH problem for Q. 

8.2.2 Construction for -1 < a < 

The above constructed Q can be used as a parametrix P for T in case a > 0. For —1 < a < 0, 
the parametrix should also have the same jump as T on {—5, 0), and we seek a 3 x 3 matrix 
valued function P such that 

1. P is analytic in Bs \ (Ai U Af U {-6, 0)) . 
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2. P has a jump P+{z) = P-{z) Jt^z) on each of the oriented contours shown in Figure [H 
right. They are given by 



hix) 



JT[X) = 

3. P{z) behaves near the origin hke: 








—x~ 


a 








+ z~ 




+ \x\ 


Qgn(A3 




as z ^ 0. (8.36) 



4. As n ^ oo, 

P(z) = Na{z){I + 0{l/n)) uniformly for z G \ (Ai U Af U (-5, 0)) , (8.37) 

where A^^^ is the parametrix built in Section [6j 

Just as in Case 1, we build P in the form ()8.22p . 

P{z) = Qiz)S{z), (8.38) 

where Q is the matrix valued function constructed by formulas (|8.33p ~ (|8.35p . and S is now 
explicitly given in each of the four components of Bs \ (Ai U Af U {—5, 0)) by 

5(z)=/-^-^z"e"(^3-A.)(.)^3^^ 

for z in the region outside the lens, (8.39) 
S{z) =1 ^^^z"e"(^3-A2){^)^ + e^'i^s-^in^) Er^^, 

for z in the upper part of the lens around Ai, (8.40) 
S(z) =1 ^^^z"e"(^=*-^2){^)^ ^_^e''{^3-^i){^)Er^^^ 

for z in the lower part of the lens around Ai. (8-41) 



Then by straightforward calculations it can again be checked that all conditions 1.-4. of 
the RH problem for P are satisfied. 

In order to unify notation for a > and — 1 < a < 0, we take as the parametrix in Bs 
in the Case 2 the matrix valued function P = QS, where 5 = / if a > 0, and S is given by 
(|839]) - (l8:in) . if -1< a < 0. 
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9 Final transformation 



We denote generically by Bs the smah disks around the branch points 0, p and q, and by P 
the local parametrices built in Bg. We define the matrix valued function R as 

\T(z)P~^(z), in the neighborhoods 
R{z) = < _ (9.1) 

\T{z)N^ {z), elsewhere. 

Then R is defined and analytic outside the real line, the lips and of the lenses and 
the circles around the three branch points. If a > 0, the jump matrices of T and coincide 
on Ai and A2 and the jump matrices of T and P coincide inside the three disks with the 
exception of the interval (0,5) in Case 1, and (—(^,0) in Case 2. It follows that R has an 
analytic continuation to the complex plane minus the contours shown in Figure [101 




Figure 10: Jump contours for the RH problem for R, when a > 0: Cases 1 (top) and 2 
(bottom). 

We find that R satisfies the following RH problem, that we describe explicitly only in the 
Case 1 (Case 2 is similar): 

1. i? is analytic outside of the contours in Figure [TOl 

2. R has a jump R+{z) = R-.{z) jRiz) on each of the oriented contours in Figure [T0| with 
jump matrix 

jR{z)=N^{z)jT{z)N-\z), Af UAfu(5,p-5)U(g + ,5,oo), (9.2) 

3r{z) = N^{z)p-\z), z G dBs, (9.3) 
jR{z)=P{z)jT{z)p-\z), zG(0,5). (9.4) 
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Figure 11: Jump contours for the RH problem for R, when — 1 < a < 0: Cases 1 (top) and 2 
(bottom). 



3. R{z) = 1 + 0{l/z) as z ^ oo. 

Note that it is only after this final transformation that the RH problem is normalized at 
infinity. Item 3. follows from ()5.3p and (16. 3p and the definition (j9.ip of R. 

If — 1 < a < 0, the situation is even simpler, since now R has an analytic continuation 
to the complex plane minus the contours shown in Figure [TTl so that only jumps (|9.2p - ()9.3p 
remain. By ()8.20p and ()8.36p . R{z) is at most O d^^l") as z — > 0, so that the singularity at 
is removable. 

From the matching conditions for the local parametrices it follows that 
jR{z) = I + C>{\/n) as n ^ oo uniformly for z on the boundary of the disks. 
If a > 0, for X in the interval (0, 5) (in Case 1) or (—5, 0) (in Case 2), we have for some c > 0, 

j^(x) = / + 0(x"e-'="). 
On the remaining contours we have for some c > 0, 

j^(z) = I + 0(e-^"l^l) asn^oo. 
We can use standard arguments (see e.g. [7]) to conclude that 

uniformly for z in the complex plane outside of these contours. Then by Cauchy's theorem 
also 
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Thus, we obtain the following estimate which will be useful in the next section 

R-\y)R{x) = I + R-\y) (i?(x) - R{y)) =I + . (9.7) 

10 Proofs of the theorems 

The proofs of Theorems I2.4fl2^ are based on the asymptotic analysis of the kernel Kn{x,y). 
If we use ()2.24p and follow the steps of the RH steepest descent analysis, we find that for 
x,y > and x,y £ Ai, 



Kn{x,y) = —j^ r (0 wi{y) W2{y)) Y;\y)Y+{x) ( 



27ri{x — 


y) 


1 




27ri(x — 


y) 


1 




2Tri{x — 


y) 


1 




27ri(x — 


y) 



l_-^(0 y"e-"^ o)^+'(2/)^+(x) 

/gnAi,+ (a;)N 

(O y"e-"^2,+(j/) o) C/+^(y)C/+(x) 





gnAi,+ (x) 

^_g-nAi,+ (y) ya^-n\2,+ {y) q) T^^(y)r+(x) I X-"e"^2,+ (x) 1 (^^q 



This will be our basic formula for the kernel. 







Proof of Theorem \2.4\ We take x and y in the interior of Ai, and we may assume that the 
circles around the branch points are such that x and y lie outside of these disks, so that 

T{x) = R{x)N^{x) , T{y) = R{y)N^{y) . 

Thus, by ([921) 

T-\y)T+{x) = N-\{x)R-\x)R+{y)N^,+ {y) 

= N-Xix) (/ + o(^))iVa,+ (y) 
= I + O {x — y) as y — > x , 

and also 

r+\y)r+(2;) x"" o\=I + 0{x-y) as y ^ x . 

.0 1/ \0 1/ 



Taking into account that on Ai both Ai and A2 are purely imaginary on Ai and A2+ = Ai+ 
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on Ai, we can rewrite (jlO.ip as 

gftj Im Ai__|_ (x) 



e 



KJx, y) = (_e-™in^^i.+(f) (.niimM^+iy) q) (J + O (x - y)) I e-™i°i^i.+W 

2'Ki[x — y) ' 







1 /'g™Im{Ai^+(j/)-Ai,+ (x)) _ g-mIm{Ai,+ (s/)-Ai,+ (a:)) _|_ ^j-^ _ y\\ 

C — v) \ / 

+ 0(1), asy^x, (10.2) 



27ri(x — y) 

sin (nIm(Ai,+ (y) - Ai,+ (x))) 



7r(x - y) 

where 0(1) holds uniformly in n. Now we let y ^ x. Using (j4.9p and the L'Hopital rule, we 
get that 

71 Tl 

Kn{x,x) = ImCi,+ (x) + 0(1) = - |ImCi,+ (x)| +0(1), n ^ oo , 

TT vr 

(see e.g. ()4.8p ). and so 

lim -Knix,x) = -|ImCi +(x)| . 
n^oo n vr 

If X € M+ \ Ai, then it can be proved analogously that 

lim —Kn{x, x) = . 

n— >oo 77, 

This proves that the limiting mean density of paths exists and is supported on This 
proves Theorem 12.41 □ 



Proof of Theorem \2.7[ Let x* S {pj^{t),q{t)), where p+{t) < q{t) are the end points of the 
interval Ai, described in Theorem l2.41 Then p{x*) > 0, where p is the density given in ()2.25p . 
For given x,y £M, we take 

* I * I y 



np{x*)^ " np{x*) 

Then for n large enough, we have Xn,yn S (j'+(*)) ^(0)) so that (|10.2p holds. Then by Taylor 
expansion, 

Im(Ai,+ (2/n) - Ai,+ (Xn)) = (yn " x„)Im Ci,+ (a^* ) + 0{yn - Xnf 

, / *^^ , 1 



T^jx-y) ^ ^ ( 1 



and therefore 



np[x*j 7r[x — y) 

= ^^^^f^ + 0(l/n), 
7T[X - y) 

which proves Theorem 12. 7[ □ 
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Proof of Theore'm \2.8l Take c = f'{q) where / is the conformal map from (jT.ip . For y G M 
we put Xn = q + ^^2/3 ^iid Un = Q + ■ This imphes that 

n2/3/(x„,) ^ X, n2/3/(yn) ^ 2/- 

If a;,y < 0, then we still can apply (|1U.1|) . but now, for n large enough, XmUn belong to the 
small disk Bs around q, so that 

T{Xn) = R{Xn)P{Xn) 

= R{xn)E{xn)-^ {n^'^fixrS) diag (^^^-"/^ef x^/^^-f (A2(x„)-Ai(x„))^ ^ 

and similarly for T{yn)- Therefore, 

and 

(_e-'^Ai,+(s/„) yag-nA2,+(y„) o)r+^(y„) 

= y;^/2e-f (-1 1 0) (n2/3/(y„)) i?"! . 

As in [71 Section 9], we can show that 

E-\yn)R~\yn)R{Xn)E{Xn) ^ /• 



Thus, 



hm -\-Kn{Xn, Vn) = ^--j^ r ("1 1 O) ^l\x)^! +{y) \ 1 

_ Ai(x)Ai^(j/) - Ai^(x)Ai(j/) 
x-y 

Similar calculations give the same result if x and/or y are positive. 

The scaling limit near p in case t < t* follows in a similar way. □ 

Proof of Theorem \2.^A Now we assume t > t* so that we are in Case 2. For x and y are in 
the (5-neighborhood of 0, we use the expression (jlO.ip for Kn{x,y) with T = RP = RQS, 
where S = / in case a > 0, or S" is given by (j8.40p in case — 1 < a < 0. In both cases it 
follows that 

^g-nAi,+ (s/) ^Qg-nA2,+ (y)^ Sl^^{y) = (e~"'^i'+ "^^z) yag-"^2, + (y)^ ^ 







/ g"Ai,+ (x) \ 








\ 1 




^ / 
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so that by (flOTI) 

(gnAi,+ {x) \ 
^-agnA2,+ (x) 1 _ (IQ 3) 

Let now x,y > be arbitrary. Let c = — /'(O) > where / is the conformal map from 
(|8-34p and take x„, = y„ = ^ so that 

nV(a;n) ^ -2;/4, nVCyn) ^ -2//4 

as n ^ oo. Then for n large enough, we have that Xn and yn are in the (5- neighborhood Bs 
of 0, so that we can use (I10.3P with x and y replaced by Xn and y„. We then have 

R+{Xn)Q+{Xn) = R{Xn)E{Xn)^+{n'^f{Xn)) 

X diae f e"''*/^x~"/^e5^'^^'+^^"^"'^^'+^^"^^ g-a'ri/2^a/2g-f (A2,+ {x„)-Ai,+ (x„)) 
and similarly for Rj^{yn)Q+[yn)- Thus, 

(gnAi,+ (a;„) \ / 

and 

(_e— yag-nA2,+(y„) q) (y„) 

= y-/2e-f (Ai,+ (y„)+A2, + (,„)) (_g-a.^/2 ^^^2 q) (n2/(2/„)) E~\yn)R^\yn) • 

Then it may be shown (see (|9.7p and [7]) that 
and we arrive at 

lim -^Knixn^Vn) = ^^-T T i-T'^ (-e-"-/2 e«W2 q) $;i(y/4)$+(x/4) ( e'^^/^ 

n-^oo cn^ 1tti[x — y) \xJ ' 
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To evaluate this further, we first note that by definition of ^, 

(^_^-ani/2 q) (y/4)$+ (x/4) | g'^^^/^ 



^^-a7ri/2 ga7ri/2^ 1-1^ (-y/4)^'„ (-x/4) ^ 



ga7ri/2 
^—aiTi/2 



^aiTi/2^ 

g— a7ri/2 

-1 



where for the last line we used the definition of ^'(C) in terms of the Hankel functions that is 
valid for — vr < arg(^ < — 27r/3. Since 

\ {h^^ + H^)) = J. 

and since the above matrices with the Hankel functions have determinant one, it follows that 
the above expression is equal to 

Using this in the expression for the scaling limit we obtain the theorem. □ 



11 Appendix: approach via equilibrium measures 

In the appendix we indicate an approach via equilibrium measures. Our starting point is 
the RH problem for X, see Proposition 13.11 Instead of the A-functions that come from the 
Riemann surface we use the so-called g-functions to make the second transformation of the 
RH problem. 

As an intermediate step we first define 

U{z) = X{z) diag(l, e2"v^/*, g-^^V^/*) (n.l) 

with the usual principal branch of the square root function. Then U satisfies the following 
RH problem. 

1. U{z) is analytic in C \ (M U ). 
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2. U{z) possesses continuous boundary values on R U A2 denoted by [/+ and [/_, and 



C/+(a;) = ?7_(x)(l + x"e Vt(i-t) * ) Ex2] 





'1 

U+{x) = U-{x) I 

,0 Ixl" 



, X E {-00, P-), 



C/+(x) = C/_(x) 







X e {p-,0), 



\ 



(11.2) 
(11.3) 

(11.4) 
(11.5) 



3. U{z) has the fohowing behavior near infinity: 
U{z) = {I + 








^ 




^1/4 

















( 






( 


Vo 













z^oo, zGC\(]RuAf), (11.6) 



4. C/(^) has the same behavior as X{z) at the origin, see (j3.20p . 

Now we consider the fohowing variational problem for two measures /xi and ^2- Minimize 



/ / loSl — - — \d^ii{x)d^i{y) - [ [ log- — - — -dfii{x)dfi2{y) 
JJ \x-y\ JJ \x-y\ 



+ 



log 



x-y\ 



dfi2ix)dfi2{y) + 



X 



2\/ax 



t{l-t) 



t 



dfiiix) (11.7) 



over all pairs {^1,112) such that 

supp(/Ui) C [0,cx)), j dm = 1, 
supp(/i2) C (-00,0], j dfi2 = 1/2, 

and 

where a is the (unbounded) measure on (— oo,0] with density 



dx vrt ' 



X £ (—00, 0]. 



(11.8) 



(11.9) 



(11.10) 
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It is possible to show that there is a unique minimizing pair (/ii,/i2)- The measures are 
absolutely continuous with respect to Lebesgue measure and their densities are related to the 
functions ^ and ^3 coming from the Riemann surface as follows 

Thus 

supp(/ii) = Ai, supp(^2) = (-00,0], supp(o- - 112) = A2, 

and the constraint pi.9p on ^2 is active only in Case 2. 

The following variational equalities and inequalities hold for certain Lagrange multipliers 
li and I2: 

2/iog|.-»i<i,,M-/iog|.-.|dM»)-j^+^{;;;; (n.12) 

This is a vector equilibrium for the pair of measures /ii and fi2, supported on and ]R__, 
respectively, with the matrix of interaction 

2 -1 
-1 2 

characteristic of a Nikishin system [13^ |46] (see [2] for a survey), but with two additional 
features: 

(i) there is an external field 

X 2JI1X 



ip{x) 



t{l-t) t 

acting on , motivated by the varying character of the orthogonality weights in (j2.20p ; 

(ii) there is an upper constraint (|11.9p originated in the fact that W2/W1 is the Cauchy 
transform of a discrete measure on ]R_, see (I3.2p . The upper constraint (Ill.lOp is equal 
to the limiting distribution of the points p.Sp that are related to the positive zeros of 
the Bessel function Jq,. 

We introduce the ^-functions 

gj{z) = j \og{z - s)dfijis), j = 1, 2, (11.14) 

and define the transformation 

U{z) = C7„diag (e""'Sl,e"'2) C/(z)diag (^e~''^'^'^'^~^'\e''^3iiz)-g2{z) ^^n{g2{z)-h)^ ^^^ -^^^ 

where li and I2 are the constants from (jll.l2p and pi.lSp and Cn is a constant matrix (see 
the first matrix in the right-hand side of (j4.18p ). Then U satisfies a RH problem. 
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1. U{z) is analytic in C \ (M U A2 )• 

2. U{z) possesses continuous boundary values on M U denoted by {/+ and and 



;7+(x) = U-{x) 








„n(gi+(x)-gi_(x)) 





'1 

C/+(x) = C/_(x) I 



, X £ (-00, 



U+{x) = U-{x) 





^0 
/l 

Q ^Ain\ax\^/'^ /t^-n(g2+{x)-g2~{x)) q 

apn{gi+{x)-g2+{x)~g2-{x)+l2) ^-Mn\ax\^/'^ /t ^n{g2+{x)-g2-{x)) 




= C/_(z) (^/ + e±"^^z-"e-^"("^)'^'/*e"(292W-9iW-'2)^23) ^ G . (11.19) 
3. [/ (z) has the following behavior asz— >oo,zGC\(MU A^): 

/l 



U{z) 



I + O 



'10 
z^l^ 
.0 



75 71^ 
VO 73^ 73, 









a/2 











(11.20) 



4. C/(^) has the same behavior as X{z) at the origin, see (j3.20p . 

Due to the equilibrium conditions we have that the jump (|11.16p simplifies on the interval 



Ai to 



' ^-■n.(gi+{x)-gi-(x)) q\ 



U+{x) = U^{x) 



Ij 



(11.21) 



A calculation that uses the fact that /i2 = cr on (p_ , 0) shows that the diagonal entries of the 
jump matrix (jll.lSp on (p-,0) are equal to 1, so that 



U+{x) = U^{x) (^I +\x\''e''^3l+{^'^-92+{x)-g2^{x)+l2) ^ ^ {p-,Q), 



(11.22) 



with an off-diagonal entry that is tending to as ?i ^ cxo. Of course the jump (jll.22p is only 
relevant in Case 2. 

We can then go on by opening a lens around Ai as discussed in the main part of the text. 
We conclude this appendix by giving the relation between the ^-functions and the A- 
f unctions coming from the Riemann surface. We have 



gi{z) = \i{z) - i 



gi{z) - 92(2) 

92{z) 



-A3(^) + 



z 



2\faz 



t{l-t) t 
z 



t{l-t) t 



(11.23) 
(11.24) 

(11.25) 
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with constants £i, £2, and £3 appearing in (j4.13p - (|4.15p . These relations and (j4.13p - (j4.15p 
show that 

..|.) = log.- ''-"<'t°"-"' +Om. (11.26) 



z \ z 



1, t + 4a(l-t) (1 - t)(t + a(l - 1)) , ^ , . „s 



as z — > 00. 
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